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ABSTRACT 


Dynamical scattering of binaries and triple systems of stars, planets, and com¬ 
pact objects may produce highly inclined triple systems subject to Kozai-Lidov (KL) 
oscillations, potentially leading to collisions, mergers, Type la supernovae, and other 
phenomena. We present the results of more than 400 million gravitational scattering 
experiments of binary-binary, triple-single, and triple-binary scattering. We compute 
the cross sections for all possible outcomes and explore their dependencies on incoming 
velocity, mass, semi-major axis, and eccentricity, including analytic fits and discussion 
of the velocity dependence. For the production of new triple systems by scattering we 
find that compact triples are preferred, with ratios of outer to inner semi-major axes 
of ~few-100, flat or quasi-thermal eccentricity distributions, and flat distributions in 
cosine of the mutual inclination. Dynamically formed triples are thus subject to strong 
KL oscillations, the “eccentric Kozai mechanism,” and non-secular effects. For single 
and binary flyby encounters with triple systems, we compute the cumulative cross sec¬ 
tion for changes to the mutual inclination, eccentricity, and semi-major axis ratio. We 
apply these results to scattering events in the field, open clusters, and globular clus¬ 
ters, and explore the implications for Type la supernovae via collisions and mergers, 
stellar collisions, and the lifetime and dynamical isolation of triple systems undergoing 
KL oscillations. An Appendix provides an analysis of the velocity dependence of the 
collision cross section in binary-single scattering. 


1 INTRODUCTION 

Gravitational scattering events are common in globular clus¬ 
ters and are also dynamically impor tant for man y 
terns in open clusters and the field (Hills fc Day 
Although two-body scattering can be studied analyti¬ 
cally, three-body scattering is to o complex to permit gen¬ 
eral, practical analytic results (IPoincar 1 Il892l ; ISundmanl 
1 1907f) . The three-body scattering problem was therefore 
not studi ed in detail until the development of comput- 


lanysys- 

iJtlizl). 


ety of astrophysical phenomena fe.e., Holman et al. 

1997|: 

Ford et al.ll2000: iMiller & Hamiltonll2002l: Blaes et al. 

20021; 

Wu & Murray 20031: Wen 20031: 

Fabrvckv & Tremaine 

2003 ; 

Wu et al. 20071: Ivanova et al 

20081: Perets & Fabrvckv 

20091: Naoz et al.l 

201 ll: ThomDsonl 1201 ll: iNaoz et al. 

20121: Katz & Dong 

20121: ShaDDee & Thompson 20131: 

Antonini et al.l 2014: 

Antognini et al. 20141). Furthermore, it 


is now known that triple systems are not rare in our Galaxy; 
demographic surveys of the field hav e revealed that triples 


ers 

fl98: 


e.g., ISaslaw et al.l Il974l ; iHeggiel 1 19751 : iHut fc Bahcalll 
Since then binary-single and binary-binary scat- 


const i tute 10% of all s t ellar s ystems (|l>ujiuennoy_&^Iay(x 
1991; iRaghavan et al. 201(]| ; iTokovininl 120141 ; ISana et al 
201411 . 


1982 


Mikkolal 19831: iHilld 1 

199ll: Valtonen & Mikkola understand in detaiHiowAriples interact with other stellar 

1991 

Sigurdsson & Phinnevl 

1993!: Bacon et al. 1999: systems. Observations of the Tarantula Nebula have demon- 

Fregeau et al. 20041. and the references therein), but lit- strated that binary interactions strongly affect the observed 


tie attention has been devoted to the scattering of triple 
systems or to the orbital characteristics of triple systems 


2002 

Ivanova et al. 2002 

: Leigh et al. 201ll: Leigh & Geller 

2012 

Moeckel & Bonnet 

120131: ILeigh & Gelled 2013. 

20151). 

Yet due to the Kozai-Lidov (KL) mechanism (Lidov 

19621: 


iKozail Il962l ) , it is now appreciated that the dynamics 
of triples may play an important role in a wide vari- 


may contribute to the observed multiplicity fraction as well. 
Moreover, the large multiplicity fraction of high mass stars 
indicates that triple d ynamics may be even more impor¬ 
tant for these systems (ILeigh fe Gelled 120131 ): for example, 
scattering of triples with binaries may be a n important for¬ 
matio n mechanism for quadruple systems. IRaghavan et al.l 
d20ldh showed that quadruple systems are nearly as com- 
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2 Antognini & Thompson 


mon as triple systems and many even higher -order systems 
have been discovered (e.g., IKoo et abll2014l b but it is un¬ 
clear what fracti on of such systems are for med in situ ver¬ 
sus dynamically (iGoodwin fc Kroupall2005lb Quadruple and 
higher-order systems may be even more dynamically impor¬ 
tant than triples since numerical and semi-analytic exper¬ 
iments have demonstrated that they exhibit stronger KL 
oscilla tions than triples over a wider region of parameter 
space dPeicha et al.ll2013l : lHamers et al.ll20liih . 

The perturbative influence of interloping stars on hi¬ 
erarchical triple systems may have an important effect on 
the long-term evolution of triples. The KL mechanism can 
drive the inner binary to very high eccentricities, but the 
maximum eccentricity reached is sensitive to the orbital pa¬ 
rameters, in particular the mutual inclination and the outer 
eccentricity. In the standard quadrupole order KL mecha¬ 
nism the tertiary can drive the inner binary to arbitrar¬ 
ily large eccentricities over a narrow range of inclinations 
near 90° (the exact value depends on the particular system). 
However, at octupole order, in the so-called “eccentric KL 
mechanism,” the range of inclinations over which the ter¬ 
tiary can drive the inner binary to extreme eccentricities is 
significantly increased, especially if the eccentricity of the 
orbit of the tertiary is large. Thus, perturbations to the in¬ 
clination or eccentricity of the tertiary may produce much 
stronger KL oscillations for star-star or star-planet systems. 
If these KL oscillations drive the inner binary to sufficiently 
high eccentricities, t he components of the inner binary will 


20031: Fabn 

mkv & Tremaind 20071: IPerets & FabrvckvH2009l: 

Naoz et al. 

201 ll. 20121: Naoz & Fabrvckv 20141'), potentially 


explaining the observ ation that nearly all close binaries 
d Tokovinine^al. 2006|) a nd cer t ain subsets of warm and hot 


Jupiters dWu fc M urray! 20031: IWu et al. 2007 : Naoz et alJ 


I201ll.l2012l : ISocrates et al.lteOlllDong et alJl20li T have ter¬ 

tiary companions. 

Tidal interactions may not alway s occur in highly in¬ 
clined triple systems; iLi et al.l d2014h showed that some 
coplanar systems can be driven to very large eccentrici¬ 
ties. Moreover, certain systems undergo changes in eccen¬ 
tricity so rapidly that the angular momentum of the inner 
orbi t can change by an order of magnitude in a single or¬ 
bit dAnt onini fc Peretsll2012l:lKatz fc Don p2012l ; ISetoll201^ : 


lAntonini et al.ll2014l : Antognini et al.ll2014l ') . These rapid ec¬ 

centricity oscillations present the possibility that the com¬ 
ponents of the inner orbit may be driven to merger more 
rapidly than the tidal circularization process can circularize 
the orbit and quench KL oscillations. The components of 
the inner binaries of such systems will then not merge gen¬ 
tly, but will instead collide head-on. These head-on collisions 
may produce a variety of unusual astrophysical phenom¬ 
ena depending on the objects comprising the inner binary. 
Most notably, if the inner binary components are two white 
dwarfs, these collisions may produce a Type la supernova 
(SN la; |ThompsorJ|201l|; |Katz_&M3ong| 2012|: lHamers et al.l 


120131 : [Kushnir et al.l 20131 : Prodan et al. 2013lb 

The possible connection between triple dynamics and 
SNe la should be explored because, despit e their crucial role 
in constraining cosmol ogical parameters dRiess et al .11 199 
IPerlmutter et ablll999l h it is unknown whether the progen¬ 
itor syste ms consist of one white dwarf (the sin gle degener¬ 
ate model; IWhelan fc Ibenlll973l : lNomotolll982l l or two (the 


doub le degenerate model; llben fc Tutukovlll984l : IWebbinl3 
Il984h . Although observati onal evidence currently favors the 
doub l e degenerate mod el dHowe H l201ll : iMaoz fc Mannuccil 
120121: iMaoz et al.l l2014lj . it is unclear how to drive white 
dwarf binaries to merge at a rate l arge enough to be consis¬ 
tent with the observed SN la rate dRuiter et al.ll2009l . l201li j. 
One approach is by driving the binary to high eccentricity 
through the perturbative influence of a tertiary via the KL 
mechanism. At high eccentricities the inner binar y wou ld 
then emit much more gravita t ional radia tion dBlaes et alJ 
120021: iMiller fc Hamiltonl 120021: IWeidl2003l) . ther eby leadi ng 
to mo re rapid coalescence of the WD-WD binary dThompsonl 
Eon|). However, it is unclear how to prevent the inner binary 
from merging or colliding while still on the main sequence. 
KL oscillations would bring the two stars into tidal con¬ 
tact and tidal circularization would then shrink the orbit, 
greatly increasing the semi-major axis ratio and “freezing” 
the in ner binary at an inclination close to the critical Kozai 
angle dFabrvckv fc Tremaindl2007l ') . These two effects would 
effectively shut off further KL oscillations and prevent the 
inner binary from merging when the stars evolve into white 
dwarfs. The dynamical formation of high-inclination triples 
from scattering may be one means to circumvent this diffi¬ 
culty. 

Scattering of high-multiplicity systems may also be 
one channel to produce free-floating planets. Microlensing 
studies have revealed that there may be as many as two 
free-floating planets for ever y bound planet in the Galaxy 
dZapatero Osorio et alJ l2000l : ISumi et abl 1201 il l . Moreover, 
there is some evidence that it is difficult for planet-planet 
scatterin g to produce free floatin g planets in the required 
numbers dVeras fc Ravmondll2012l l . Dynamical scattering of 
field stars off of binary stars may be an important compo¬ 
nent to the rate of planet ejection. 

Scattering events involving high-order stellar systems 
could also be a source of stellar collisions. Stellar collisions 
are estimated to occur at a rate of ~0.5 yr _1 in the Galaxy 
dKochanek et al.ll2014l l and it is unknown whether these col¬ 
lisions are due to KL oscillations, scattering events, or a 
combination of the two. Even when scattering does not lead 
to a stellar collision the scattering event will change the or¬ 
bital parameters of the system. This will lead to evolution 
in the distribution of orbital parameters of the triples in the 
population. This evolution would be strongest in globular 
clusters due to their high stellar densities and large ages. 

This paper has three goals. The first is to perform ideal¬ 
ized binary-binary, triple-single, and triple-binary scattering 
experiments in order to derive general relationships between 
the cross sections of various outcomes with the initial orbital 
parameters and the incomin g veloc ity. I n this way our p aper 
is similar to iHut fc Bahcalll dl983l l and iMikkola d 19831 1 ex¬ 
cept that we study higher-order scattering. The second goal 
of our paper is to determine the orbital parameters of triple 
systems after scattering events. The final goal of our paper 
is to apply the cross sections we derive in several contexts. 
In particular we estimate the rate of formation of WD-WD 
binaries in highly inclined triples, the ejection rate of plan¬ 
ets due to scattering events, the stellar collision rate, and 
the lifetime of high inclination triples. To meet these goals 
we have performed over 400 million scattering experiments. 

Details of the numerical methods of this paper are pre¬ 
sented in Section [2] We study triple scattering in detail in 
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Dynamical formation & scattering of triples 3 


Section[3]by performing numerical experiments and compar¬ 
ing them with analytic approximations in Section[4] We then 
discuss the distribution of orbital parameters of dynamically 
formed triples in Section [5] We discuss the implications of 
these results for SNe la, the longevity of triple systems un¬ 
dergoing KL oscillations, and provide estimates for the rate 
of planetary ejection and stellar collisions in Section [6] We 
conclude in Section [7] 


2 NUMERICAL METHODS 

We use the open so urce Fewbody suite to perform our scat¬ 
tering experiments dFregeau et al.l[2004l j. Fewbody is opti¬ 
mized to numerically compute the dynamics of systems with 
small numbers of components (N < 10). Fewbody uses the 
ordinary differential equations library of the GNU Scientifi c 
Library (GSL) for its underlying integrator (iGou 53 ill). 
The GSL ordinary differential equations library supports six 
integration algorithms with adaptive time steps, from which 
we use eighth-order Runge-Kutta Prince-Dormand integra¬ 
tion. We find, however, that the choice of integration algo¬ 
rithm makes little difference to the results of the calculations 
because GSL’s adaptive time steps are chosen to target a 
specified relative and absolute accuracy (10 -14 in our ex¬ 
periments) regardless of the algorithm used. Fewbody also 
supp orts the use of Kustaanheimo -Stiefel (KS) regulariza¬ 
tion dKustaanheimo fe StiefeJ Il965l i. a coordinate transfor¬ 
mation which removes the singularities in the gravitational 
force present in ordinary IV-body integration. Throughout 
this paper we use KS regularization as it has the particular 
advantage of making eccentric orbits much easier to com¬ 
pute. 

2.1 Notation 

Throughout this paper we use the same notation as 
iFregeau et alJ d2004l 'l and Fewbody in which each star in 
a system with n stars is labelled with a unique index run¬ 
ning from 0 to n — 1. Bound pairs are denoted by square 
brackets surrounding the pair and collisions between two 
stars are denoted by colons between the pair. For example, 
a hierarchical triple with an unbound interloping star is no¬ 
tated [ [0 1] 2] 3 (stars 0 and 1 form the inner binary, star 
2 is the tertiary, and star 3 is the interloping star), and a 
hierarchical triple in which the two stars of the inner binary 
have collided is notated [0:1 2]. 

Orbital parameters (e.g., semi-major axis or eccentric¬ 
ity) of the component binaries of the system are notated by 
subscripts in a top-down fashion. Leftmost indices in sub¬ 
scripts represent the outermost binaries of the separate hi¬ 
erarchies, and rightward indices represent inner binaries in 
those hierarchies. For example, in a triple-binary scattering 
event, the semi-major axis of the outer binary of the triple 
is ai, the semi-major axis of the inner binary of the triple is 
an, and the semi-major axis of the interloping binary is a 2 . 
Similarly, the mass of the second star in the innermost bi¬ 
nary of the triple is mm, the mass of the tertiary is 77112 , and 
the mass of the first star of the interloping binary is 77121 . 
We also combine masses in the subscripts. Thus the total 
mass in the inner binary of the triple is mn = mm + mm 
and the total mass of the triple is mi = mn + 77112 . 


For clarity we occasionally refer to the orbital param¬ 
eters of hierarchical triples using the subscripts “in” and 
“out.” Thus the eccentricity of the inner orbit of a triple 
may be referred to equivalently by en or ei n and the eccen¬ 
tricity of the outer orbit by ei or e ou t- 

2.2 Cross sections 

Throughout this paper we provide cross sections for the out¬ 
comes of scattering events involving triple or liigher-order 
hierarchical systems. Our treatment of these cross section s 
and their uncertainties follows that of iHut fe Bahcall (11983h . 
We use the usual definition of the cross section for a partic¬ 
ular outcome, X: 

,2 Tlx /1 \ 

U A' — 7T0 max , (ij 

Titot 

where & max is the maximum impact parameter in a set of ex¬ 
periments, Tlx is the number of experiments with outcome 
X, and ntot is the total number of experiments performed. 
Throughout this paper we present our results using normal¬ 
ized cross sections, a. The cross sections from triple-single 
scattering experiments are normalized to the area of the 
outer binary of the triple: 

<7 = (triple-single). (2) 

7ra, 

The cross sections for binary-binary experiments are nor¬ 
malized to the sum of the areas of the two binaries and 
similarly in triple-binary experiments the cross sections are 
normalized to the sum of the areas of the outer binary of 
the triple and the interloping binary: 

<7 = . f 7 —(binary-binary, triple-binary). (3) 

TT(ai + ai) 

There are two uncertainties in the calculation of the 
cross sections. The first is the statistical uncertainty due to 
the finite number of experiments performed^ 

(J y 

AstatOW = - (4) 

\Jn X 

The second source of uncertainty is due to the fact that cer¬ 
tain systems require prohibitively long computation times 
to resolve. For example, while it is impossible for a sin¬ 
gle star scattering off of a triple system to produce a sta¬ 
ble quadruple system JChazvlll929l ; lLittlewoodlll952l : iHeggiel 
119751 : iHeggie fe Hutll2003h . certain systems can enter into a 
“metastable” state where a quadruple system is produced 
that takes an exceedingly long time to dissociate. Such un¬ 
resolved systems produce a separate systematic uncertainty 
in the calculation of the cross sections given by 

A ,2 7l U nres / r , 

sys&X — 7rO ma x 5 (o) 

Tltot 

1 In cases where the fraction of outcomes with the outcome X is 
close to unity or zero it is better to use a confidence interval like 
the Wilson score interval rather than the normal approximation 
in order to ca pture the asymmetry of the statistical uncertainty 
(|Wilsonl Il927f ) . However, we perform enough experiments that 
the difference between the Wilson score interval and the normal 
approximation is almost always negligible. Where it is not, we use 
the Wilson score interval. 
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4 Antognini & Thompson 


where n unres is the number of experiments with unresolved 
outcomes. 

It is important to note that the systematic uncertainty 
is completely asymmetric and only serves to increase the 
cross section of any outcome. That is to say, our estimates of 
the cross sections are lower bounds on the true cross sections. 


2.3 Initial conditions and halting criteria 


Each scattering experiment begins with the interloping sys¬ 
tem at a large, but finite, distance from the target system. 
We choose the initial separation between the two systems to 
be the distance at which the tidal force on the outer binary 
of the target system is some small fraction, <S, of the relative 
force between the two components of the outer binary when 
at apocenter. Hence, the initial separation, r, is given by 


Ttid 

E re i 


= 5, 


( 6 ) 


where 


Ftn = 


2G(m 11 + mi 2 )m 2 

p 3 


a(l + e), 


(7) 


and 


F le i = 


Gmnmi2 

[a(l + e)f 


( 8 ) 


where a and e refer to the semi-major axis and eccentric¬ 
ity, respectively, of the outer binary of the target system. 
T hroughout this paper we use $ = 10 -5 , the same choice 
as lFregeau et al.l (I2004l l. Smaller choices of 5 do not change 
the results but increase the running time. The interloping 
system is then analytically brought along a hyperbolic orbit 
from infinity with the 6 max and velocity at infinity, Voo , that 
has been fixed for that particular experiment. 

It is also essential to choose the appropriate b m ax when 
computing cross sections. Too small a choice of 6 max will re¬ 
sult in in an underestimate of the cross section, whereas too 
large a choice of i> max will result in few experiments produc¬ 
ing outcomes of interest, thereby leading to large statistical 
uncertainties. We ado pt the choice of 6 max similar to that of 
I Hut fe BahcaH (Il983ll of: 


b 


max 


( 4Ucrit 

\ U in f 



m, 


(9) 


where Vi n f is the incoming velocity at infinity and u cr it is the 
critical velocity at which the total energy of the system is 
zero. For binary-binary scattering, u cr it is 


l^crit, bin-bin 


G(m i +m 2 ) (m u m 12 m 2 im 22 ^ 


mim 2 


ai 


a 2 


( 10 ) 


for triple-single scattering, v cr it is 


.,2 _ G(aiminmn 2 + anmnmi 2 )(mi + m 2 ) 

^crit, trip-sing 


aianmim2 
and for triple-binary scattering v C rit is 


(ii) 


2 


^crit, trip-bin 


G(mi + m 2 ) 

CL 1 CL 11 CL 2 vn 1 77i 2 


x [aiia 2 mnmi 2 + ai(a 2 minmn 2 + onm 2 im 22 )]. 


( 12 ) 


We will frequently refer to velocities in terms of the normal¬ 
ized velocity, v, which is the velocity, v, scaled to Vcrit, 


^crit 


(13) 


We test our choice of & max by computing the cross sections 
for four categories of outcome in triple-single scattering over 
a broad range of fe max . We demonstrate that the cross sec¬ 
tions are well converged in Fig. [T| 

The calculation halts when one of four conditions is 
met: (1) if the system consists of some number of stable 
hierarchic al systems as determined b y the Mardling stability 
criterion (iMardling fe Aarsethll200lll : 


— (1-ei) >2.8 
an 


1 _|_ 77111 ^ 1 + e i 

m 12 ) yjl - ei 


2/5 


1 0.3. 

1- 1 

7T 


(14) 

(where the mutual inclination, i, is in radians), and all the 
systems are far enough apart that the ratios between their 
tidal forces to their relative forces is less than 5; (2) the sys¬ 
tem integrates for 10® times the initial orbital period of the 
outer binary of the target system; (3) the system integrates 
for one hour; or (4) the total energy or angular momentum 
of the system changes from its initial value by more than 
one part in 10 3 . If conditions (2) or (3) are met the outcome 
is classified as unresolved and included in our systematic 
uncertainty^ 


3 SCATTERING EXPERIMENTS 

Hierarchical triples can be formed through scattering in one 
of three ways: (1) binary-binary scattering, (2) scattering of 
triples, or (3) scattering of liigher-order systems. Binary¬ 
binary scattering has been studied extensively, but only 
rarely have any studies addressed the formation of triple 
systems and, moreover, scattering of higher-order systems in 
general has received little attention at all. We study triple 
and binary scattering in detail in this section. We do not 
study the scattering of higher-order systems (e.g., triple- 
triple, quadruple-single) be cause such scattering events are 
rare in most environments (ILeigh fe Gellerll2013l l. 

We compute the cross sections for the outcomes of both 
single and binary stars scattering off of hierarchical triples 
along with binaries scattering off of binaries in the point 
mass limit using Newtonian gravity. The masses of all stars 
in the system are equal unless we explicitly vary the mass 
of one star of the inner binary of the triple. In an individual 
scattering event, the semi-major axis and eccentricity of all 
component binaries (both binaries of the triple and, in the 
case of triple-binary scattering, that of the interloping bi¬ 
nary) in the system are fixed, but the argument of pericenter 


2 IPortegics Zwart fe Boekholtl |j2014ll have argued that typical 
energy conservation standards in IV-body studies are extremely 
conservative and that maintaining energy conservation to better 
than one part in 10 is sufficient to preserve the statistical proper- 
tie s of dynamical experiments. These r esults have been confirmed 
bv lBockholt fe Portegies Zwart) 120 1 5| ). Since very few of our ex¬ 
periments violate our energy conservation standard, this finding 
does not change our results, but does imply that we may have 
slightly overestimated our systematic uncertainties. 
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Figure 1. Convergence of the scattering cross sections with 6 max for triple-single scattering. We performed 3 x 10 5 scattering experiments 
at each choice of 6 max . Left panel: All masses are equal point masses and the orbits of the initial triple were circular with a semi-major axis 
ratio of 10. The cross sections for four categories of outcomes are presented: exchanges (black circles), single ionization (blue triangles), 
double binary formation (red squares), and exchange-ionization (orange stars). (See Section 13.21 for a description of the categories.) Right 
panel: Same as the left panel, but the mass of one star in the inner bin ary is set to be a test particle. The choice for 6 max we adopt in 
this paper (dashed line) is given by equation (9} and is similar to that of lHut &; Bahcalll d 1983ll . Smaller Umax’s do not probe the full area 
over which non-flyby outcomes occur, whereas the statistical power of larger Umax’s is reduced because of the small fraction of non-flyby 
outcomes. The choice we adopt is converged but is small enough to have strong statistics. 


and the mean anomaly are chosen from a uniform distribu¬ 
tion between 0 and 2n. The component binaries are then 
oriented randomly by pointing each binary’s angular mo¬ 
mentum vector toward a randomly chosen point on a sphere. 
We present a sample binary-binary, triple-single, and triple¬ 
binary scattering experiment projected in the x?/-plane in 

Fig.H 

We first perform scattering experiments on a model sys¬ 
tem in Section T3. II and then vary several of the initial orbital 
and physical parameters of this system one at a time to de¬ 
termine the dependence of the cross sections on these pa¬ 
rameters in Section m 

3.1 Cross sections of model systems 

Our model system for triple-single scattering is a moder¬ 
ately hierarchical triple system with semi-major axis ratio, 
a = ai/on = 1 0. This choice of a is larg e enough that the 
triple is stable (IMardling fc Aarsethll2001l ). but is still small 
enough that there will be dynamical interactions between 
the inner and outer binaries. In this model system, the inter¬ 
loping star has an incoming velocity of Voo/ucrit = 1, which 
is approximately the velocity dispersion of stars in the thin 
disk of the Milky Way (~40 km s _1 ) relative to the critical 
velocity of a triple system consisting of three 1 Mr stars 
with ai = 10 AU and an = 1 AU feinnev fc Merrifieldl 
119981, p. 656). We take both orbits to be circular. 

Our model system for triple-binary scattering is identi¬ 
cal, but with the interloping star replaced by an interloping 
binary with semi-major axis equal to the outer semi-major 


axis of the triple, ai = a\ = 10an. The incoming velocity is 
again taken to be u C rit, but with v cr it now calculated using 
equation m instead of equation CD- Again, all stars are of 
equal mass, so now the mass ratio of the interloping system 
to the triple is 2/5 instead of 1/3. 

Lastly, our model system for binary-binary scattering 
is identical to the case of triple-binary scattering but with 
the inner binary of the triple replaced by a point mass of 
1 Mg. The incoming velocity is taken to be v cr it, but with 
Ucrit calculated using equation (flOll. In this case the mass 
ratio of the interloping system to the target system is 1. 

We perform 10 6 scattering experiments for each model 
system. We present the cross sections for all possible out¬ 
comes with their statistical and total uncertainties in Ta¬ 
bic [T] for the triple-single case, Table [5] for the triple-binary 
case, and Table [3] for the binary-binary case. 

Because there are a large number of possible outcomes 
in binary-binary, triple-single, and triple-binary scattering 
(22, 22, and 161, respectively), many of which are qualita¬ 
tively similar, we group these outcomes into broad classes. 
In the case of triple-single scattering we define six classes: 
(1) flybys, in which the hierarchical structure of the sys¬ 
tem remains the same, although the orbital parameters may 
have changed; (2) exchanges, in which the interloping star 
replaces one of the stars in the triple; (3) double binary for¬ 
mation, in which one star of the triple is ionized and binds to 
the interloping star; (4) single ionization, in which one star 
from the triple is ionized, leaving a binary and two unbound 
stars; (5) full ionization, in which all stars become unbound 
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Figure 2. A sample binary-binary, triple-single, and triple-binary scattering experiment projected onto the xy- plane. The incoming 
velocity is 0.1 v cr n in the case of triple scattering and 0.2t) cr ;t in the case of binary-binary scattering. The semi-major axis ratio of the 
triple in the case of triple scattering is 10 and the semi-major axis of the incoming binary is equal to the semi-major axis of the outer 
binary in the case of triple-binary scattering. In the case of binary-binary scattering the semi-major axes are equal. In the case of triple 
scattering the triple approaches from the left and the interloping systems approaches from the right. Left panel: In this binary-binary 
scattering event one star from each binary combine to form a new binary, and the two remaining stars are ejected. Middle panel: In this 
triple-single scattering event one star from the inner binary of the triple is ejected and the interloping star is captured as a tertiary. Right 
panel: In this triple-binary scattering event one star from the inner binary is ejected, leaving two binaries and a single star. 


Table 1. Normalized cross sections for the outcomes of triple¬ 
single scattering (see equation[2]for the normalization). The initial 
conditions are described in Section [3] and the initial hierarchy is 
[ [0 1] 2] 3. We present both the statistical uncertainty and the 
systematic uncertainty. Note that the systematic uncertainty only 
represents an uncertainty toward larger cross sections. That is, 
the cross sections presented are lower limits and may be larger 
by the systematic uncertainty—see Section 12.21 We only present 
cross sections for which the statistical uncertainty is less than 50 
per cent. The cross section for the outcome equal to the initial 
hierarchy (i.e., a flyby) is not well defined and so is not included. 
Such interactions can change the orbital parameters of the triple 
system, however, and are discussed in Section 15.11 The outcome 
classes are defined in Section 13.21 We also note if the scattering 
event produces a new triple with a hierarchy distinct from the 
original hierarchy. 


Triple-single 


Outcome 

a 

Astat<5" 

Outcome class 

[0 1] 3 2 

1.309 

0.013 

Single ionization 

[0 2] 1 3 

0.092 

0.003 

Single ionization 

0 [1 2] 3 

0.080 

0.003 

Single ionization 

0 [1 3] 2 

0.076 

0.003 

Exchange ionization 

[0 3] 1 2 

0.074 

0.003 

Exchange ionization 

[[0 1] 3] 2 

0.033 

0.002 

Exchange, new triple 

[[0 3] 2] 1 

0.027 

0.002 

Exchange, new triple 

0 [[13] 2] 

0.021 

0.002 

Exchange, new triple 

[0 1] [2 3] 

0.017 

0.001 

Double binary 

[0 2] [1 3] 

0.011 

0.001 

Double binary 

3 [[12] 0] 

0.010 

0.001 

Scramble, new triple 

[0 3] [1 2] 

0.010 

0.001 

Double binary 

[[0 2] 1] 3 

0.008 

le-03 

Scramble, new triple 

0 1 [2 3] 

0.005 

8e-04 

Exchange ionization 

AsysO" 

0.032 

0.002 



from each other; and (6) a scramble, in which the tertiary 
exchanges with one of the stars of the inner binary. 

In the case of binary-binary scattering we define six 
classes: (1) flybys, defined as in the triple-single case; (2) 
exchanges, in which one star from each binary exchanges 
places with the other; (3) triple formation, in which a stable 
triple is formed, leaving a single unbound star; (4) single 
ionization, in which one star is ionized, leaving a binary 
and two unbound stars; (5) full ionization, defined as in the 
triple-single case; and (6) exchange + ionization, in which 
one star from each binary binds to form a new binary, and 
the two other stars remain unbound. 

In the case of triple-binary scattering we define seven 
classes: (1) flybys, defined as in the triple-single case; (2) 
exchanges, in which one star from the interloping binary ex¬ 
changes with one star from the interloping binary; (3) dou¬ 
ble binary formation, in which one star from the triple is 
ionized, resulting in two binaries and an unbound star; (4) 
binary disruption, in which the two stars of the interloping 
binary become unbound; (5) triple disruption, in which the 
three stars of the triple become unbound from each other; (6) 
quadruple formation, in which one star from the interloping 
binary becomes bound to the triple, forming a quadruple; 
(7) and a scramble, defined as in the triple-single case. 

Note that in triple-single and triple-binary scattering 
several of these outcome classes produce a ‘new’ triple (i.e., 
the final hierarchy differs from the original hierarchy). (In 
binary-binary scattering triple formation is its own class 
and is mutually exclusive with the other classes.) The cross 
sections for these classes for triple-single, triple-binary, and 
binary-binary scattering are displayed in Table[4] We include 
there the cross section for new triple formation, though for 
triple-single and triple-binary scattering it is not indepen¬ 
dent of the cross sections for other outcome classes. 

In the case of triple scattering the systematic uncer¬ 
tainty is dominated by marginally unstable triples (accord¬ 
ing to Equation 11-11) for which the integration time exceeded 
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Table 2. Normalized cross sections for the outcomes of triple¬ 
binary scattering (see equation [3] for the normalization). The ini¬ 
tial conditions are described in Section l3l and the initial hierarchy 
is [ [0 1] 2] [3 4]. We present both the statistical uncertainty 
and the systematic uncertainty. Note that the systematic uncer¬ 
tainty only represents an uncertainty toward larger cross sections. 
That is, the cross sections presented are lower limits and may be 
larger by the systematic uncertainty—see Section l2.2l We only 
present cross sections for which the statistical uncertainty is less 
than 50 per cent. Flyby outcomes are not included. Such inter¬ 
actions can change the orbital parameters of the triple system, 
however, and are discussed in Section 1 5.1 1 The outcome classes 
are defined in Section |3.2| We also note if the scattering event 
produces a new triple with a hierarchy distinct from the original 
hierarchy. 


Triple-binary 


Outcome 

<7 

Astatd 

Outcome class 

[0 1] 4 2 3 

1.110 

0.016 

Double ionization 

[0 1] [3 4] 2 

0.990 

0.015 

Double binary 

[[0 1] 2] 4 3 

0.911 

0.015 

Binary disruption 

[0 1] 3 [2 4] 

0.072 

0.004 

Double binary 

[0 1] 4 [2 3] 

0.070 

0.004 

Double binary 

[0 2] 14 3 

0.048 

0.003 

Double ionization 

0 [1 2] [3 4] 

0.048 

0.003 

Double binary 

0 [1 2] 4 3 

0.048 

0.003 

Double ionization 

[0 2] 1 [3 4] 

0.045 

0.003 

Double binary 

[[0 1] 4] 3 2 

0.025 

0.002 

Bin. disruption, new triple 

0 [1 3] 2 4 

0.025 

0.002 

Double ionization 

[0 4] 12 3 

0.022 

0.002 

Double ionization 

[[0 1] 3] 4 2 

0.021 

0.002 

Bin. disruption, new triple 

[0 3] 12 4 

0.021 

0.002 

Double ionization 

0 [1 4] 2 3 

0.018 

0.002 

Double ionization 

[[0 1] 3] [2 4] 

0.011 

0.002 

Exchange, new triple 

[[0 1] 4] [2 3] 

0.011 

0.002 

Exchange, new triple 

[[0 2] 1] [3 4] 

0.006 

0.001 

Scramble, new triple 

[[0 3] 2] 14 

0.005 

0.001 

Bin. disruption, new triple 

[0 4] [1 3] 2 

0.005 

0.001 

Double binary 

[3 4] [[1 2] 0] 

0.005 

0.001 

Exchange, new triple 

0 1 [2 3] 4 

0.004 

le-03 

Double ionization 

0 1 2 [3 4] 

0.004 

le-03 

Double ionization 

[0 3] [1 2] 4 

0.004 

le-03 

Double binary 

[[0 2] 1] 3 4 

0.004 

9e-04 

Bin. disruption, new triple 

[[0 4] 2] 13 

0.004 

9e-04 

Bin. disruption, new triple 

3 [[1 2] 0] 4 

0.004 

9e-04 

Bin. disruption, new triple 

0 [[1 3] 2] 4 

0.003 

9e-04 

Bin. disruption, new triple 

[0 2] [1 3] 4 

0.003 

9e-04 

Double binary 

[0 4] [1 2] 3 

0.003 

9e-04 

Double binary 

[0 3] [1 4] 2 

0.003 

9e-04 

Double binary 

[[0 4] 3] 12 

0.003 

8e-04 

Bin. disruption, new triple 

0 1 [2 4] 3 

0.003 

8e-04 

Double ionization 

0 [[1 4] 2] 3 

0.003 

8e-04 

Bin. disruption, new triple 

[0 2] [1 4] 3 

0.002 

7e-04 

Double binary 

0 [[1 4] 3] 2 

0.002 

7e-04 

Bin. disruption, new triple 

[0 3] 1 [2 4] 

0.002 

7e-04 

Double binary 

[0 2] [[1 4] 3] 

0.002 

7e-04 

Exchange, new triple 

[0 4] 1 [2 3] 

0.002 

7e-04 

Double binary 

[0 4] [[1 3] 2] 

0.001 

6e-04 

Exchange, new triple 

[[0 3] 2] [1 4] 

0.001 

6e-04 

Exchange, new triple 

0 [1 3] [2 4] 

0.001 

6e-04 

Double binary 

[[0 4] 2] [1 3] 

0.001 

5e-04 

Exchange, new triple 

[0 3] [[1 4] 2] 

0.001 

5e-04 

Exchange, new triple 

[[0 2] 3] 14 

0.001 

5e-04 

Bin. disruption, new triple 

0 [[1 3] 4] 2 

0.001 

5e-04 

Bin. disruption, new triple 

3 [[1 4] 0] 2 

0.001 

5e-04 

Bin. disruption, new triple 

[[0 3] 4] 12 

0.001 

5e-04 

Bin. disruption, new triple 

[[0 4] 1] 3 2 

0.001 

5e-04 

Bin. disruption, new triple 


Table 3. Normalized cross sections for the outcomes of binary¬ 
binary scattering. The initial conditions are described in Section 
[3]and the initial hierarchy is [0 1] [2 3]. The cross sections have 
been normalized to the sum of the areas of the two binaries. We 
present both the statistical uncertainty and the systematic uncer¬ 
tainty. Note that the systematic uncertainty only represents an 
uncertainty toward larger cross sections. That is, the cross sec¬ 
tions presented are lower limits and may be larger by the system¬ 
atic uncertainty—see Section 12.21 We only present cross sections 
for which the statistical uncertainty is less than 50%. The cross 
section for the outcome equal to the initial hierarchy (i.e., a flyby) 
is not well defined and so is not included. 


Binary-binary 


Outcome 

(T 

Astatd 

Outcome class 

0 1 [2 3] 

1.335 

0.016 

Single ionization 

[0 1] 3 2 

1.310 

0.016 

Single ionization 

[0 2] 13 

0.829 

0.013 

Exchange ionization 

[0 3] 12 

0.823 

0.013 

Exchange ionization 

0 [1 3] 2 

0.807 

0.013 

Exchange ionization 

0 [1 2] 3 

0.800 

0.013 

Exchange ionization 

[0 2] [1 3] 

0.152 

0.005 

Exchange 

[0 3] [1 2] 

0.150 

0.005 

Exchange 

[[0 3] 2] 1 

0.003 

8e-04 

Triple formation 

[[0 1] 2] 3 

0.002 

6e-04 

Triple formation 

[[0 3] 1] 2 

0.002 

6e-04 

Triple formation 

0 [[1 2] 3] 

0.002 

6e-04 

Triple formation 

0 [[1 3] 2] 

0.002 

6e-04 

Triple formation 

[[2 3] 0] 1 

0.002 

6e-04 

Triple formation 

[[0 2] 1] 3 

0.002 

6e-04 

Triple formation 

2 [[1 3] 0] 

0.002 

6e-04 

Triple formation 

[[0 2] 3] 1 

0.001 

5e-04 

Triple formation 

[[0 1] 3] 2 

0.001 

5e-04 

Triple formation 

3 [[1 2] 0] 

0.001 

5e-04 

Triple formation 

0 [[2 3] 1] 

0.001 

5e-04 

Triple formation 

Asysd 

0.002 

0.001 



the one-hour CPU time limit; 0.0263 per cent of the triple¬ 
single scattering experiments failed to complete, and of these 
81 per cent failed to complete due to the CPU time limit 
(the rest failed because they violated the energy conserva¬ 
tion limit). Because the Mardling stability criterio n of equa¬ 
tion m is not a hard boundary (|Petrovicbll2015h it is pos¬ 
sible that these triples are, in fact, stable and should be 
classified as flybys. If, however, they are unstable on longer 
timescales than our integrations permit they should be clas¬ 
sified as single ionization events since no other outcome is 
possible. Thus the systematic uncertainty of triple scatter¬ 
ing should be considered as an uncertain contribution to 
the cross section for single ionization. In the case of binary¬ 
binary scattering the systematic uncertainty is dominated 
by systems that violated the energy conservation criterion 
during an extremely close passage. This is a much rarer oc¬ 
currence than the formation of a marginally unstable triple 
during a triple scattering event, so the systematic uncer¬ 
tainty of binary-binary scattering is much lower than that of 
triple scattering (0.0012 per cent of systems failed to resolve, 
of which all were due to violations of the energy conservation 
limit). 
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Table 4. Cross sections for outcome classes in triple-single, triple¬ 
binary, and binary-binary scattering at v = 1. Cross sections for 
“Exchange + ionization” are not included in either “Exchange” or 
“Single ionization.” However, the cross sections for new triples are 
not independent of the other cross sections. (I.e., they are a sum 
of subsets from the other classes.) See Section 13.21 for definitions 
of the outcome classes. 


Triple-single 


Outcome class 

<7 

Astatd 

Exchange 

0.082 

0.003 

Single ionization 

1.481 

0.013 

Full ionization 

0.000 

0.000 

Double binary formation 

0.038 

0.002 

Exchange + ionization 

0.155 

0.004 

Scramble 

0.018 

0.001 

New triple 

0.100 

0.003 

Asys<7 

0.032 

0.002 

Triple-binary 

Exchange 

0.036 

0.003 

Quadruple formation 

0.000 

2e-04 

Scramble 

0.006 

0.001 

Binary disruption 

0.996 

0.016 

Double binary formation 

1.250 

0.017 

Triple disruption 

1.303 

0.018 

Full ionization 

0.000 

2e-04 

New triple 

0.127 

0.006 

Asys<5" 

0.045 

0.003 

Binary-binary 

Exchange 

0.302 

0.008 

Triple formation 

0.020 

0.002 

Single ionization 

2.646 

0.023 

Full ionization 

0.000 

0.000 

Exchange + ionization 

3.259 

0.025 

Asys<7 

0.002 

0.001 


3.2 Dependence on initial parameters 

We next explore the dependence of the cross sections on the 
initial parameters of the system. In particular we separately 
vary the semi-major axis ratio, the incoming velocity, the 
eccentricities, and the masses. 

3.2.1 Semi-major axis ratio 

We first vary the semi-major axis ratio from a = 02/01 = 
10 0,6 — 10 2 . At choices of a much below our lower limit the 
triples become unstable and at choices of a much larger 
than our upper limit the computational time becomes pro¬ 
hibitively long because the time step is set by the orbital 
period of the inner binary, but the crossing time is set by 
the size of the outer orbit. We hold all other parameters (e.g., 
eccentricities, masses) fixed to their values in the model sys¬ 
tem (Section l3.ll) . The normalized cross sections for binary¬ 
binary, triple-single, and triple-binary scattering are shown 
in Fig. [3] Single ionization is the dominant outcome of triple¬ 


single scattering for all a. At the low-a end this is because 
the triples are only marginally stable so minor perturba¬ 
tions from the interloper tend to ionize one member of the 
triple. At the high-a end this is because the tertiary be¬ 
comes weakly bound to the inner binary. Weak perturba¬ 
tions from the interloper are therefore likely to ionize the 
tertiary. Other classes of outcomes generally require some 
interaction with the stars of the inner binary. The probabil¬ 
ity of such an interaction scales with the orbital area of the 
inner binary, so we observe that a oc a -2 . 

For the case of triple-binary scattering we fix the initial 
value of ai to be equal to 02 as we vary the semi-major axis, 
a. At large a, binary disruption is the dominant outcome 
whereas at small a double binary formation dominates. In 
the low-a limit the triple comes closer to instability, so single 
ionization from the triple is also likely to occur in addition 
to the disruption of the interloping binary, leading to double 
binary formation. In the large a case the scattering problem 
approaches that of binary-binary scattering in which one bi¬ 
nary is more ma ssive than the other . Earlier experiments 
have shown (e.g.. iFregeau et al.ll2004l ) that the less massive 
binary is more likely to be disrupted, in accordance with 
our results here. In this limit we might naively expect no 
dependence on a since we fix ai = 02 . Instead we observe 
a dependence of <7 oc a -1 . This is because as a increases 
the incoming velocity (which, in the large-a limit is com¬ 
parable to the orbital velocity of the inner binary of the 
triple) becomes larger relative to the orbital velocities of the 
outer binary and interlopin g binary. In the limit of large in¬ 
coming velocities a oc v~ 2 (IHut fe Bahcalllll983l) . and since 
v oc a _1 ' /2 , we have that <7 oc an oc a -1 . For the same 
reason the scaling of the single ionization cross section in 
triple-single scattering exhibits the same dependence. 

3.2.2 Eccentricity 

We next vary the eccentricity of binaries in the systems, 
both individually and together while holding constant all 
other parameters (e.g., semi-major axis ratio and masses). 
To increase the range over which we can vary the outer ec¬ 
centricity in triple scattering we set the outer semi-major 
axes to 20 AU (i.e., we use a of 20 instead of 10 as before). 
We vary the outer eccentricity of the triples (ei) from 0 to 
0.7 and the eccentricity of the inner binary (en) from 0 to 
0.98. Even with the larger value of a at outer eccentricities 
larger than ~0.7 the triple violates the Mardling stability 
criterion (equation |T4j) . In the case of triple-binary scatter¬ 
ing we fix ei = e 2 . In the case of binary-binary scattering we 
vary the eccentricities of either one binary or both binaries 
from 0 to 0.98. The cross sections for triple-single, triple¬ 
binary, and binary-binary scattering are shown in Fig. [I] 
In the case of binary-binary scattering the cross sections 
are completely independent of the eccentricity. In the case 
of triple scattering the cross sections are independent of the 
inner eccentricity and are generally independent of the outer 
eccentricity, although there is a modest increase in the cross 
section for single ionization in the case of triple-scattering 
and in the cross section for double binary formation in the 
case of triple-binary scattering. 

The physical mechanism behind these dependences is 
twofold. First, there is some baseline cross section for 
the outcome which is independent of any eccentricity. 
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Binary-binary scattering Triple-single scattering Triple-binary scattering 
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Figure 3. Cross sections of various classes of outcomes as a function of semi-major axis ratio, a , for binary-binary (left panel), triple¬ 
single, (middle panel), and triple-binary (right panel) scattering. In the limit of large semi-major axis ratios for triple-single scattering, 
interactions involving the inner binary should scale like the area of the inner orbit (i.e., an inverse square dependence on the semi-major 
axis). Indeed, the cross sections for scrambles and exchange + ionization (both of which require the interloper interact with the inner 
binary) are consistent with the inverse square dependence shown (black dashed line). Single ionization is the dominant outcome of 
triple-single scattering at all semi-major axes. At small a the triple becomes unstable to ionization from weak perturbations due to the 
interloper and at large a the tertiary is weakly bound relative to the inner binary so weak perturbations from the interloper lead to 
ionization. See Section 13.2.11 for a discussion of the origin of the a dependence of different cross sections. 


I Hut &; Bahcalll (jl983l ) argue that this is because the cross 
sections for a particular orbit are only dependent on the 
average velocity of the orbit, which is independent of the 
eccentricity as a result of the virial theorem. This baseline 
cross section applies to all outcomes. Second, in the case of 
triple scattering there is an additional contribution to the 
cross section of ionization-like outcomes (i.e., single ioniza¬ 
tion in the case of triple-single scattering and double binary 
formation in the case of triple-binary scattering) due to per¬ 
turbations to the eccentricity of the outer orbit of the triple, 
ei. If ei exceeds some critical eccentricity, e C rit, the triple 
will violate the Mardling stability criterion and will become 
unstable. As we show in Section GEU and in particular in 
panels b) of Fig. 1101 the logarithm of the cross section for 
cumulative changes in the eccentricity follows the inverse of 
the Gompertz function: 

Aei = exp (—ci exp (02 lncf)) = exp (—cicr C2 ), (15) 

where ci and C 2 are positive constants. This implies that the 
cross section to undergo an ionization-like reaction takes the 
form 

a = do + (-uf 1 ln(e cr it - ei)) 1/c2 , (16) 

where do is the baseline cross section described above. Tak¬ 
ing 6 C rit from the Mardling stability criterion (e cr it ~ 0.733), 
we find an excellent match between this functional form and 
the cross sections for ionization-like outcomes in Fig. [I] Best¬ 
fitting parameters are shown in Table [5] 

Because ei n does not influence the stability of the triple, 
we do not see a similar effect when ei n is varied. This result 
that the cross sections are independent of the eccentricities, 
except when related to the stability of the triple, is consis¬ 
tent with the results of binary-single scattering experiments, 
which have also fou nd that cross sections are independent 
of eccentricity (e.g., iHut fe Bahcalll fl983T ). We note, how¬ 
ever, that is only the cross sections of strong interactions 


Table 5 . Best-fitting parameters of the eccentricity cross sections 
in Fig.[4]to the inverse Gompertz function. See equation 1161 1 for 
the definition of the parameters. 


Scattering type 

panel 

<5-0 

Cl 

C2 

Triple-single 

b 

0.35 

6.21 

0.57 

Triple-single 

c 

0.34 

6.05 

0.60 

Triple-binary 

b 

0.23 

3.45 

0.35 

Triple-binary 

c 

0.14 

4.88 

0.39 


that are independent of eccentricity. iHegaie fc Rasiol (i l996i f 
found that the magnitude of small secular perturbations to 
the eccentricity of a binary from a distant flyby is propor¬ 
tional to e\J 1 — e 2 . 


3.2.3 Mass ratio 

We show in Fig. [5] the effect of changing the mass of one 
component of the system in binary-binary, triple-single, and 
triple-binary scattering. In the case of triple scattering the 
mass of one component of the inner binary of the triple 
is varied. All other parameters (e.g., a, eccentricities, and 
v ) are held fixed. The mass of the one star is varied from 
0.03 Mq to 40 Mq. For small mass ratios (i.e., as the test 
particle limit is approached) the cross section for all out¬ 
comes in triple scattering increases by roughly an order of 
magnitude. In the high mass limit single ionization becomes 
the dominant outcome of binary-binary scattering, rising 
roughly with the square root of mass, and binary disrup¬ 
tion becomes the dominant outcome of triple-binary scat¬ 
tering with roughly the same mass dependence. Conversely, 
although single ionization is the dominant outcome at high 
mass for triple-single scattering, the cross section decreases, 
nearly with the square of the mass. 
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Figure 4. Cross sections of outcome classes for triple-single (top row), triple-binary (middle row), and binary-binary (bottom row) 
scattering as a function of eccentricity. For triple scattering we vary the inner eccentricity alone (panels a), the outer eccentricity alone 
(panels b), and the inner and outer eccentricities simultaneously (panels c). In the case of triple-binary scattering we set the eccentricity 
of the incoming binary to be equal to that of the tertiary. In the case of binary-binary scattering we vary the eccentricity of one binary 
only (left panel) and both binaries simultaneously (right panel). In binary-binary scattering the cross sections are independent of the 
eccentricities. In triple scattering the cross sections are independent of the inner eccentricity and are independent of the outer eccentricity 
for small outer eccentricities. At large outer eccentricities there is a small increase in the cross section for single ionization in triple-single 
scattering and double binary formation in triple-binary scattering. This is because these triples start closer to the Mardling instability 
boundary, so small perturbations to the outer eccentricity can destabilize them. We fit these cross sections with the functional form of 
equation 1161 ) (dashed lines) and find excellent agreement. The best-fitting parameters are shown in Table [5] 
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iHutl (Il983f) derived the mass dependence of the ioniza¬ 
tion cross section for binary-single scattering in the large 
v limit and found it to be proportional to m~ 2 when the 
most massive body resides in the binary and proportional 
to m 2 when the most massive body is the interloping star 
(see equation □Hi. In the large mass limit of binary-binary 
and triple-binary scattering, the mass dependence should be 
similar to the case of binary-single scattering with a massive 
interloping star. Likewise, the case of triple-single scattering 
should exhibit a similar mass dependence to binary-single 
scattering with a massive binary. However, we are not in 
the high velocity limit (v = 1) so we expect the mass depen¬ 
dencies to soften, which is what we observe. Nevertheless, 
the trend of increased ionization cross sections for low-mass 
interloping binaries and decreased ionization cross sections 
for low-mass interloping single stars holds even in the mod¬ 
erate velocity regime. 


3.2.4 Incoming velocity 

We finally examine the effect of the velocity at infinity of 
the interloping system. We measure these velocities relative 
to the critical velocity of the system (see equations [Tl] and 
O and explore a range of two dex centered around u C rit- 
It is difficult to explore a much larger range because veloci¬ 
ties much larger than 10w C rit yield very small cross sections 
and therefore have large statistical uncertainties. Velocities 
much smaller than 0.1v cr it require an excessive amount of 
computing time because the low total energy of the system 
means that it takes a long time for the system to random 
walk to a region of phase space where one object has a suffi¬ 
cient amount of energy to escape and leave behind a stable 

system. _ _ 

The cross sections are presented in Fig. 161 a (1 198311 
found two high-velocity limits in binary-single scattering: 
the cross section for ionization, which has a dependence 
of <7 oc v~ 2 and the cross section for exchange, which has 
a dependence of o oc v~ 6 . Although there are many more 
possible outcomes in binary-binary, triple-single, and triple¬ 
binary scattering, the behavior of the cross sections is qual¬ 
itatively similar to binary-single scattering, in that at high 
velocity the outcomes may be grouped into those with a v~ 2 
dependence (“ionization-like” outcomes), and those with a 
steeper dependence (“exchange-like” outcomes). In some 
cases exchange-like outcomes follow the same v ~ 6 depen¬ 
dence as in binary-single scattering (e.g., triple formation in 
binary-binary scattering and quadruple formation in triple¬ 
binary scattering), but in others the velocity dependence is 
shallower (e.g., exchanges in binary-binary scattering, which 
display a v ~ 4 dependence). We discuss these results in more 
depth in Section EM 


4 ANALYTIC APPROXIMATIONS 

We here develop analytic approximations to the cross sec¬ 
tions computed in Section [3] We first review binary-single 
scattering, and then generalize the theory of binary-single 
scattering to four-body scattering. 


4.1 Three-body scattering 


The analytic theory of binary-single scat tering of p oint 
masses was developed bv lHeggia dl975f) and lHutl ( 1983l l and 
has been found to agree with numer ical scattering experi¬ 
ments quite well dHut fe Bahc.alllFl983l '). Although the deriva¬ 
tion is complicated, the results are straightforward to sum¬ 
marize. In binary-single scattering of point masses only three 
classes of outcomes are possible: (1) a flyby, (2) an exchange, 
or (3) an ionization. It is not possible to form a stable triple 
from binary- s ingle scattering (fc hazvdli 9291 : lLittlewoodlll952l : 
iHeggiel Il975l : IHeggie fc Huff 12003 L The c ross section of a 
flyby is not well defined, but iHutd dl983h explicitly calcu¬ 
lated the cross section for the other two classes in the high 
velocity limit (v 1) and found them to be 

40 ml 1 


3 miimufmn + m .12 + m 2 ) v 2 ’ 


v > 1 (17) 


and 


20 m 2 (m + m 2) 4 1 
3 m2(2m + m 2 ) 3 v 6 


v > 1 


(18) 


where the l atter is only valid for equal mass binaries, m = 
mn = mi 2 . IlliiTl dl983ll found that both cross sections are 
independent of eccentricity. 

In the low velocity limit (f) 1) the cross sections can¬ 

not be computed with precise numerical coefficients. Never¬ 
theless, the velocity scaling in the low velocity case is simple 
because ionization is not possible and the only mechanism 
to change the cross section is gravitational focusing. This is 
because at low velocities the speed of the interloping star 
when it is close enough to interact strongly with the binary 
will be very close to the escape speed. Thus for any low 
velocity system, the incoming velocities when the system 
begins strong interactions will be nearly the same, and the 
cross section will scale as 

a oc i<l. (19) 

v 

The cross sections may be combined by taking the re¬ 
ciprocal of the sum of the reciprocals of the cross sections in 
the two extreme limits. This is because different physical ef¬ 
fects limit the rate of production of particular outcomes and 
the combined effect of rate limiting processes is generally the 
reciprocal of the sum of the reciprocals. Because the mass 
terms are in general not known, we introduce normalization 
factors (<5"a , etc.), which can be determined empirically. 
The cross section for exchange is then 



^ex & a &b 


The validity of this functional form is demonstrated in Fig.Q 
where we fit equation I® to numerical exchange cross sec¬ 
tions (blue dashed line) for binary-single scattering of equal 
mass stars. 

The cross section for full ionization is necessarily zero 
for v < Ucrit ■ For velocities slightly in excess of v cr it (v = 
u C rit + v2e, where e <C 1), the cro ss section for full ionization 
in equal mass systems scales as (IHeggie fe Sweatmanlll99lh 

*.«! OC£ VT3 - 1 Oc(f)-l) ( ^- 1)/2 . (21) 


Therefore the general cross section for ionization is 
1 (D _ 1 )(l-^3)/2 ^ 

^ion &c O'd 


( 22 ) 
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Binary-binary scattering Triple-single scattering Triple-binary scattering 
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Figure 5. Cross sections of outcome classes for binary-binary (left panel), triple-single (middle panel) and triple-binary (right panel) 
scattering as a function of the mass ratio of one component of the binary (in the case of binary-binary scattering) or of one component 
of the inner binary (in the case of triple scattering) to all other components of the system. The masses of all other components of the 
system remain fixed and equal to each other. In the high mass limit the more massive system tends to disrupt the incoming system. This 
leads to an increase in the cross section for single ionization in binary-binary and triple-binary scattering. In the case of triple-single 
scattering all cross sections decrease because the incoming system cannot disrupt. 



10 


2 Triple-binary scattering 


10 

10 ° 

10" 1 

10~ 2 

10“ 3 


Tl. 


l h 


o 

t T i ♦ • •• * I . *i 


0 exchange 
^ quadruple formation 

# scramble 

A binary disruption 

★ DB formation 
double ionization 
full ionization 

O new triple 


f' 1 


ft 




'i 


* I- - A 

l u*4 

■ •.« * 


1 0 -4L_ 

10" 1 


10 u 


10 1 


Figure 6. Cross sections of outcome classes for binary-binary (left panel), triple-single (middle panel) and triple-binary (right panel) 
scattering as a function of the incoming velocity. The velocities are written relative to the critical velocity of the system (see equations 1111 
and 1121 . Outcomes which require an ionization from one system follow a v ~ 2 dependence at high velocities, whereas those requiring an 
exchange follow a steeper dependence. The kink in the cross section for single ionization is explained in Section 14.21 


Note, however, that the exponent on the tr c term is only 
valid in the case of approximately equal masses. If one object 
contains more than ~t /2 of the mass, the dominant orbital 
configuration leading to full ionization will change, leading 
to a different exponent dSweatmanll2007| ) . ISweatmanl (120071 ) 
lists exponents for other configurations, but this exponent 
could also be determined empirically. The validity of the 
combined cross section for ionization is also illustrated in 
Fig. [7] (red dashed line). 

4.2 Four-body scattering 

Four-body scattering can occur through binary-binary or 
triple-single scattering^ Because the number of objects in- 

3 In principle a simultaneous interaction between a binary and 
two independent single stars or a simultaneous interaction be- 


teracting is the same in both cases, they generally may be 
studied in the same framework. Although there are 22 dis¬ 
tinct possible outcomes of four-body scattering, we here ex¬ 
amine the velocity dependences of three broad classes of out¬ 
comes: (1) ionization, including single and full ionization; (2) 
exchange; and (3) triple formation. 

It is instructive to study the velocity dependence of ion¬ 
ization in the extreme semi-major axis ratio limit. If the 
semi-major axes between the two binaries are very different 
or the semi-major axis ratio of the triple is very large, then 
the problem can be considered to be a three-body scattering 
problem in two limits: (1) at low velocities the smaller binary 
can be considered a point mass and (2) at high velocities one 
component of the larger binary can be ignored. There are 

tween four single stars are also possible, but such events are rare 
and so we do not study them here. 
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Figure 7. Normalized cross sections for exchange (blue circles) 
and ionization (red squares) in equal-mass binar y-single scatter¬ 
ing a s a function of incoming velocity (cf. Fig. 5 of Hut &; Bahcalj 
ll983T) . The dashed lines are an empirical fit to equations on and 
(f22i The analytic approximation is based on a combination of 
the low velocity limits (» « 1 for exchange and i —1 < 1 for 
ionization) and the high velocity limits (v 3> 1). Values for the 
normalization parameters cr a , ©, etc. are estimated from the nu¬ 
merical data with best-fitting values of a a ~ 9.76, Ob ss 28.9, 
<t c ~ 12.8, and rr,/ ps 10.3. The good agreement between the data 
and the approximation for all v implies that the two limits capture 
all the relevant physics of binary-single scattering. 


Figure 8. Cross sections for outcomes of equal mass binary¬ 
binary scattering with a semi-major axis ratio of 100 and cir¬ 
cular orbits as a function of incoming velocity. We fit the cross 
sections for single ionization to equation 12-111 (dashed line). The 
best-fitting parameters are o a ~ 2.16 X 10 “, dp, ~ 3.47, and 
<r c ~ 8.08 X 10 2 . The ionization cross section consists of two 
components: a component similar to the full ionization cross sec¬ 
tion in binary-single scattering (curved dot-dashed line), and a 
component due to hardening of the inner binary (straight dot- 
dashed line). When summed, these two components produce a 
kink at velocities intermediate between v CT n i and v) f:r j, 2 (dotted 
lines). 


therefore two critical velocities in the four-body scattering 
problem: the critical velocity of the wide binary scattering 
off of interloping system, fi C rit,i, and the critical velocity of 
the interloping system scattering off of one component of 
the wide binary, v cr it, 2 - 

In the limit that the semi-major axis of the smaller bi¬ 
nary approaches zero, the cross section for single ionization 
will appear like the cross section for ionization in three- 
body scattering given by equation (1221) with the only dif¬ 
ference being that one object in the three-body scattering 
event has a mass twice as great as the others because it is 
a binary. The effect of this is to change the exponent on 
th e (ii — 1 ) term from \/l3 — 1 to (2 + \/T0)/8 (section 4.3 
of lHeggie fe Sweatmanlll99ll ). and to change the base from 
(4) — 1) to (v— ?)crit,i). In the limit that the semi-major axis of 
one binary approaches zero the cross section for single ion¬ 
ization approaches zero for v < v cr it,i- But the cross section 
will always be non-zero because it is possible for single ion¬ 
ization to occur by hardening the smaller binary. The cross 
section for a wide encounter with a hard binary to harden 
the bi nary enough to ionize the interloping object takes the 
form (lHeggidll975l . equation 5.44) 

dp<Tj cr jt ] i = i (23) 

where o a is a normalization factor and sub-linear terms (i.e., 
terms with only a logarithmic dependence on the velocity) 
have been neglected. The total cross section for single ion¬ 


ization in either triple-single or binary-binary scattering is 
then given by 


C^ion 


V 2 ( Ob Oc 


(24) 


In Fig. [8] we present the cross section for single ionization 
for binary-binary scattering with a semi-major axis ratio of 
100, equal masses, and circular orbits. We fit the numerical 
data to the functional form of equation (HHll and find an 
excellent match. (See the caption of Fig.[8]for the best-fitting 
parameters.) 

To study exchange and triple formation we return to 
more moderate semi-major axis ratios: unity for binary¬ 
binary scattering, and 10 for tripl e-single scattering. As dis¬ 
cussed in Section 14.11 iHutl (1 1983l l derived the high-velocity 
exchange cross section to be proportional to v~ e for binary- 
single scattering. We find from our numerical experiments in 
Section 13.2.41 that the high-velocity exchange cross section 
for exchange is shallower, being instead proportional to v~ 4 
for binary-binary and triple-single scattering (Fig. [6] black 
points). 

In binary-binary scattering new triple formation occurs 
via an exchange. We find empirically that the velocity de¬ 
pendence for the new triple cross section is proportional to 
fP 6 in the high velocity limit for binary-binary scattering 
(blue points, Fig. (9j. In triple-single scattering, exchanges 
are the dominant channel for new triple formation at low 
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velocities and velocities slightly larger than the critical ve¬ 
locity. In the high velocity limit the cross section for new 
triple formation via exchange is likewise proportional to {U 6 
for triple-single scattering. 

In the case of triple-single scattering new triple forma¬ 
tion may also occur via a scramble. In the high velocity limit 
the cross section for scrambles is proportional to v~ 2 (red 
points, Fig. El). This is because a scramble is, in a sense, an 
incomplete ionization. The interloping star imparts enough 
energy to one star of the inner binary such that its semi¬ 
major axis becomes much larger than that of the tertiary of 
the original system, but not enough energy to ionize it. At 
low velocities a scramble is a rare outcome, but due to the 
shallower velocity dependence, at sufficiently high velocities 
it becomes the dominant mechanism by which to form new 
triples in triple-single scattering. 

We may combine these velocity dependences to obtain 
the general velocity dependence of the cross section for triple 
formation. For binary-binary scattering, the cross section 
has an identical form to the exchange cross section in binary- 
single scattering, except that the break occurs near v cr it,i 
rather than -0 = 1 (this is reflected in the smaller best-fitting 
value for cp,): 



0"new trip. &a &b 


(25) 


For triple-single scattering the cross section is similar, except 
that an additional high velocity v~ 2 term must be added to 
account for scrambles: 


t± + (t+CY' 

v 2 \a b o c ) 


(26) 


We present a fit of the new triple cross section to equa¬ 
tions (HU) and (1261) in Figure El The match between equa¬ 
tion (|25l> and the binary-binary scattering cross sections is 
excellent. The match between equation (1261) and the triple¬ 
single scattering cross sections is good, but we note that at 
high velocities the scramble cross section appears to be bet¬ 
ter fit with a IT 3 / 2 dependence. We present this alternative 
velocity dependence in Figure El in gray. At present we do 
not have an analytic understanding of this shallower slope. 

While the normalization parameters a a , <76, etc., should 
be determined by htting equations 01), OH). and HU, 
for systems not too dissimilar from those presented in Fig¬ 
ures E] and El the ratios between the normalization param¬ 
eters should be approximately constant. The overall scale 
can then be estimated from the dependences of the cross 
sections in Figures HE] For systems far from the systems 
presented in Figures [5] and El the normalization parameter 
may be quite different from these estimates because they 
may depend on each other in a nonlinear way. 


5 ORBITAL PARAMETER DISTRIBUTIONS 
AFTER SCATTERING 

Scattering events can affect the orbital parameter distribu¬ 
tion of triples in two ways: by perturbing the orbital ele¬ 
ments of an existing triple in a flyby, and by creating a new 
triple system. We are particularly interested in determin¬ 
ing the distribution of orbital parameters that govern the 
strength and timescale of KL oscillations, namely cosi, a, 



Voo 


Figure 9. The velocity dependence of the cross section for new 
triple formation for binary-binary and triple-single scattering of 
our model system (see Section mi l. We fit the velocity depen¬ 
dence to equations 1251 ) and 1261 (black dashed lines). At high 
velocities the triple formation cross section from binary-binary 
scattering is proportional to 6 . At velocities slightly above the 
critical velocity, the new triple formation cross section from triple¬ 
single scattering carries the same v ~ 6 dependence, but at higher 
velocities, new triple formation is dominated by scrambles, which 
are ionization-like and would therefore carry a v~ 2 dependence. 
We note, however, that the triple-single data are better fit by a 
j)—3/2 dependence (gray dashed line). For binary-binary scatter¬ 
ing our best fitting parameters are <j a 0.567 and © rs 0.0279. 
For triple-single scattering our best fitting parameters assuming 
a v~ 2 dependence for scrambles is <j a ~ 0.247, © rs 0.0402, and 
dc Ri 0.0133. Assuming a v ~ 3 / 2 dependence for scrambles the best 
fitting parameters are o a ~ 0.257, © rs 0.044, a c ~ 8.89 X 10 -3 . 


and e ou t- We first examine the case of flybys in Section [5.II 
and then examine newly formed triples in Section [5.21 


5.1 Changes to the orbital parameters from flybys 

Cross sections for flybys are not well defined because they 
diverge for infinitesimally small perturbations. As such, we 
instead calculate the cumulative cross section, a(X > x), 
defined to be the cross section for a change in parameter X 
by at least x. From our experiments involving the model sys¬ 
tem (see Section T3.II) we select those which result in a flyby. 
We then calculate the cumulative cross section for generat¬ 
ing changes in the orbital parameters of a given magnitude 
or larger. These cross sections are shown in Fig 1101 The 
normalized cross section for inducing a change in the semi¬ 
major axis of order unity is a ss 0.54, which is very close to 
the cross section for single ionization of 1.48. That the two 
cross sections should be comparable is to be expected since 
there is only a small difference in energy in changing a by 
order unity and in ionizing the object. The logarithm of the 
cumulative cross section for changes to ei is well described 
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by the inverse of a Gompertz functiorQ (cf. equation 11511 : 

Aei = exp (—ci exp (02 In a)) = exp (—cid C2 ). (27) 

We present a fit to the cumulative cross sections in Fig. 1101 
but the difference between the two is typically smaller than 
the width of the lines (the best-fitting parameters are pre¬ 
sented in the caption to Fig. m- The cross section for large 
changes to the eccentricity is small (it C 1 for 1 — Aei <S 1) 
because to remain stable such systems need to gain a large 
amount of energy and are therefore more likely to be ionized 
than remain bound. The cross section for moderate changes 
to the eccentricity (Aei ~ 0.5) is of order unity (a ~ 1), sim¬ 
ilar to the case for a change in a. The cumulative cross sec¬ 
tion for changes to e ou t drops more steeply for Ae ou t > 0.7. 
This is due to the fact that the initial triple is unstable for 
e ou t ~ 0.7 (see eauationfldll. To produce Ae > 0.7 two things 
are required: (1) the semi-major axis ratio must increase by 
enough that at the final e ou t the triple is no longer unstable, 
and (2) e ou t must change by Ae ou t- The addition of require¬ 
ment (1) induces the sharper drop in the cross section for 
Ae ou t > 0.7. To check if the form of equation (H7]> is inde¬ 
pendent on the assumption of equal masses we repeated the 
calculations used to produce panel b) of the upper row of 
Figurc[4]and panel b) of the upper row of Figure [10] using the 
following sets of masses: (mm,mn 2 ,mi 2 ,m 2 ) £ {(4, 3, 2, 
1), (1, 2, 3, 4), and (1, 10 —6 , 1, 1)}. We found equation (l27l) 
to be an excellent fit in all cases. While this does not prove 
that equation (1271) is necessarily valid for all possible mass 
ratios, it demonstrates that the form of equation l|27|l is not 
strongly dependent on the assumption of equal masses. 

The change in the inclination is particularly interesting 
due to its relevance for inducing KL oscillations in a triple 
system. However there is a confounding effect in the cumu¬ 
lative cross section for this orbital parameter. Since some 
triples will begin at high inclinations, their inclination can 
change in isolation due to KL oscillations. This is a diffi¬ 
cult effect to disentangle from changes that are solely due 
to scattering or perturbations from interloping stars. To es¬ 
timate the magnitude of this effect, we evolve a sample of 
triple systems identical to those used in our model systems 
but without any interloping stars for the same length of 
time as our flyby calculations. We then calculate the cumu¬ 
lative distribution of the change in inclination to compare 
its magnitude to the distribution when an interloping star 
is present. Small changes in inclination can be accounted for 
almost entirely by KL oscillations. Because these scattering 
calculations complete in ~10% of the KL timescale, tKL, 
KL oscillations cannot produce large changes in inclination. 
Thus, the cumulative cross sections for large changes in in¬ 
clination are entirely due to the interloping star. In more 
compact systems, however, the KL timescale may be com¬ 
parable to the timescale of the scattering event, in which 
case either scattering or KL oscillations would be able to 
produce large changes in inclination. The distributions of 
positive and negative changes in cosi are equal to within 
the statistical uncertainty of our results. 

4 A Gompertz function is defined to be any function of the form 
f(t) = ae~ be ~ Ct , 

where a is the asymptote and b and c are positive constants. 


We additionally show the relationship between changes 
to the mutual inclination and changes to e ou t in Figure fill 
To demonstrate the magnitude of the effect of KL oscilla¬ 
tions we also include the relationship between changes to the 
inclination and e ou t from the isolated triples. Flyby events 
for which A cos i < 10~ 2 do not exhibit a strong correlation 
between the change to the inclination and e ou t. In these 
systems the dynamics of the isolated triple therefore over¬ 
whelm the contributions of the interloping star. Flyby events 
for which Acosi > 10~ 2 do exhibit a strong correlation be¬ 
tween the change in the inclination and e ou t, albeit with 
substantial scatter. This result is consis tent with other dy¬ 
namical studies (e.g., iLi fe Adamdl2015h . 

5.1.1 A note on convergence 

The issue of whether our calculations are converged in the 
maximum impact parameter, fe max , is subtler in flyby events 
than the convergence analysis presented in Section 12.31 so 
we revisit it here. In Section E3] we considered the conver¬ 
gence of the cross sections of outcomes from “strong” inter¬ 
actions (i.e., interactions for which the hierarchy changed). 
In flyby events, however, arbitrarily distant encounters can 
have arbitrarily weak dynamical effects, so our calculations 
can never be truly converged in b m ax . We can, however, de¬ 
termine the magnitude of the change of orbital parameters 
down to which our calculations are converged by calculat¬ 
ing the cumulative cross sections using a variety of 6 max ’s. 
We present the results of this calculation for changes to the 
semi-major axis ratio, a, in Figure [171 

Our results are converged down to relative changes in 
a of ~3 xlO -3 . The cross section for smaller changes to 
a diverges, however. This is due to the fact that the semi¬ 
major axes of the inner and outer orbits of an isolated triple 
will exhibit small variations over time. We demonstrate the 
scale of this variability by integrating an isolated triple in 
Figure fl3l For the triples we are studying here the scale of 
this variability is Aa/ao ~ 3 x 10~ 3 . This scale corresponds 
with the point at which the cross sections diverge in Fig¬ 
ure m Note, however, that the scale of this variability is 
dependent on the on the orbital parameters of the triples 
being studied. In particular, as q -> 00 , the relative size of 
this variability will approach zero. Calculation of the cross 
sections of changes to the orbital parameters from flybys of 
triples with very large a may therefore require a larger 6 max 
than given by equation ©. 

5.2 Distribution of orbital parameters in 
dynamically formed triples 

The scattering experiments presented in Section [3] demon¬ 
strated that some fraction of binary-binary, triple-single, 
and triple-binary scattering events produce new triple sys¬ 
tems. Because this fraction is relatively small, we here run 
three larger suites of scattering experiments to obtain bet¬ 
ter statistics on the distribution of orbital parameters of 
dynamically formed triples. One suite consists of binary¬ 
binary scattering experiments, another of triple-single ex¬ 
periments, and the last of triple-binary experiments. In the 
case of triple-binary scattering the interloping binary has 
a semi-major axis equal to the outer binary of the triple. 
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Figure 10. The cumulative cross section for changes to the orbital parameters from triple-single flybys (top row) and triple-binary flybys 
(bottom row) with incoming velocities of v = 1. We show here changes to the semi-major axis ratio (panels a), the outer eccentricity 
(panels b), and the mutual inclination (panels c). The cross sections presented are the cross sections for a change in the orbital parameter 
of that magnitude or greater. The 1 -o confidence interval is shown by the gray shaded region, although typically these uncertainties are 
smaller than the width of the lines. We fit the Aei cumulative cross sections to a Gompertz function (panels b, dotted line) and find 
an extremely close match. Except for very small or large Aei the difference between the numerical data and the fit is smaller than the 
width of the line. The best-fitting parameters are ci ~ 0.0626 and C 2 ~ 0.866 for triple-single flybys and ci ~ 0.469 and C 2 ~ 0.810 for 
triple-binary flybys. Since isolated triples can undergo changes to their inclination due to KL oscillations (an effect which is necessarily 
present when we perform our scattering experiments), we compare the magnitude of this effect to the cross sections due to scattering 
by calculating the evolution of an identical set of triples for an identical length of time (panels c, dashed line). The difference between 
the scattering cross section and the cross section from the isolated triples is shown by the thin black line with 1 -a uncertainties. KL 
oscillations are unable to produce large changes in the inclination because the flyby timescale is much shorter than the KL timescale. 


We set the semi-major axis ratio of the triples to be 10 and 
set ai = <22 in the case of binary-binary scattering. We fur¬ 
thermore set all masses to be equal and set an incoming 
velocity of v = v cr it • We do not include the orbital param¬ 
eters of triples resulting from flybys in these distributions 
as these were studied in Section 1541 To ensure that the re¬ 
sulting triples are stable we run them in isolation for 100 
outer orbital periods and reject those which do not main¬ 
tain the same hierarchy. (E.g. we reject systems which be¬ 
gin in the configuration [[0 1] 2] but end in the configu- 
rat ion [ [0 2 ] 1]). This is sim ilar to the stability test used 
bv lMardling fc Aarsethl dl999l b 

We must be careful when calculating the inclination dis¬ 
tribution. Because the calculations can often run for times 
comparable to £kl, KL oscillations can bias the inclination 
distribution towards the Kozai angles of 39° and 141°. To 
correct for this we run the triples in isolation for half the 
length of time that the scattering event was computed; this 
is approximately the length of time that the triple has ex¬ 


isted for. We then use the inclination closest to 90° in the 
inclination distribution because the triple spends most of its 
time at this inclination. (For example, an equal mass triple 
with an initial inclination of 85° spends over 80 per cent of 
its time within 10° of its initial inclination.) 

The distribution of orbital parameters is shown in 
Fig. 1141 The most important feature of these distributions 
is that scattering produces triples which are extremely com¬ 
pact, in that the ratio of the distance at periapsis of the 
outer orbit, r per i, out, to < 2 i n is quite small, typically ~10. 
Indeed, the vast majority of these triples are on the verge of 
instability. 

Because these triples are so compact, they have ex¬ 
tremely large e 0 ct- This parameter is a measure of the 
strength of octupole-order contri butions to the KL os¬ 
cillations and i s defined to be dLithwick fc Nao3 1 20 111 : 
iNaoz et al. ]») 


_ 1 
Coct = 

a 


6out 

1 P 2 
1 c out 


(28) 
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Figure 11. The relationship between changes to the mutual incli¬ 
nation and changes to e ou t from flyby events. The darkness of the 
hexagonal bins corresponds to the number of events, binned loga¬ 
rithmically. To demonstrate the magnitude of KL oscillations we 
include this relationship for several isolated triples (green points). 
For events in which there is a change to cos i less than ~10 —2 there 
is not a strong correlation between A cos i and e ou t • Changes to 
the orbital parameters are due primarily to KL oscillations rather 
than the dynamical effect of the interloping star. For events in 
which there is a change to cos i greater than ~10~ 2 there is a 
strong correlation between A cos i and Ae ou t, though with sub¬ 
stantial scatter. 


Triples with larger e oc t h ave stronger eccent ric KL oscilla¬ 
tions. In particular, iLithwick fe Naol (l201ll ) note that for 
£oct > 10“ 2 the parameter space for orbital flips (and hence 
arbitrarily large eccentricities in the inner binary) widens 
dramatically. Our results indicate that nearly all dynam¬ 
ically formed triples have e oc t > 10” 2 . We note that the 
particular triples formed from these experiments do not un¬ 
dergo typical eccentric KL oscillations because they are of 
equal mass, s o the octupole or der terms of the Hamiltonian 
vanish fe.g.. IFbrd et al.1 12000 ). Nevertheless, in extremely 
compact triples non-secular effects can dri ve the inner bi¬ 
nary to arbitrarily large eccentricities (e.g., Antonini et alJ 
120141 : iBode fc Wegdl20l3 : lAntognini et al.ll2014h . Moreover, 
as we show in Section [El these results for the dynamical 
formation of triples are generic and apply to unequal mass 
triples. 

The inclination distributions appear to be approxi¬ 
mately uniform, but are clearly biased towards inclinations 
less than 90°, i.e., prograde orbits. This is curious be¬ 
cause we began with a distribution of inclinations uniform 
in cost (i.e., as many retrograde orbits as prograde), but 
many studies have found that retrograde orbits are gener¬ 
ally more stable than prog r ade orbits (e. g.,|Harrington| [T972l : 
iMardling fe Aarsethl l200ll ; iMorais fe Giupponel 20121 ). It is 
possible that scattering is biased towards producing triples 
in prograde orbits and there remains an excess of prograde 
orbits after the unstable triples have dissociated. 

The inner eccentricity distributions have two compo- 


Figure 12. The cumulative cross section for changes to the semi¬ 
major axis ratio, a, for a variety of choices of 6 max . For changes 
larger than Aa/ao ~3x 1 0 3 our calculations using ft max = 70 
AU are converged. For changes to a smaller than ~3 X10 —3 the 
cross sections diverge with larger 6 max because a in an isolated 
triple will vary on this scale over time (see Figure fT3l) . Choices 
of &max less than 70 AU are not converged down to the limit 
of Aa/ao ~ 3 X 10 3 , resulting in a smaller asymptotic cross 
section. 
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Figure 13. The semi-major axis ratio, a over time for an isolated 
triple with initial orbital parameters of: an = 1 AU, a\ = 10 AU, 
en = e± = 0, i = 70°. The relative variation of a is ~3 xl0~ 3 , 
in agreement with the point of divergence seen in Figure fl2l 
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Figure 14. Orbital parameters of dynamically formed triples from binary-binary (solid lines), triple-single (dashed lines), and triple¬ 
binary (dotted lines) scattering. The systems consist of 1 Mq stars in circular orbits with a semi-major axis of 1 AU in the case of 
binary-binary scattering, and = 1 AU and a 0 ut = 10 AU in the case of triple scattering. We have excluded triples whose final 
hierarchy is identical to their initial hierarchy. Upper left: the inclination distribution for dynamically formed triples. We correct the 
inclinations for KL oscillations; this panel shows the minimum cos i that is achieved by KL oscillations. The distribution is approximately 
uniform in cosi. Upper middle: the distribution of inner eccentricities. Both triple-single and triple-binary scattering produce a peak 
at low eccentricities due to our choice of initial conditions of circular inner orbits. Both also produce a flat distribution of ei n whereas 
binary-binary scattering produces a distribution of ei n which is closer to thermal. Upper right : The distribution of outer eccentricities. 
Triples produced from binary-binary scattering are biased towards small e 0 ut, but triple scattering produces triples that are slightly 
biased towards large e 0 ut- Lower left: the distribution of semi-major axis ratios. Lower middle: the distribution of distances at periapsis 
relative to the inner semi-major axis. Dynamically formed triples, particularly those formed from binary-binary scattering, are extremely 
compact. Most are close to the stability boundary. Lower right: the distribution of e 0 ct, a measure of the strength of the octupole order 
term in the Hamiltonian. Generally, if e oc t > 10 —2 the EKM will operate and produce flips for a substantial range of initial inclinations. 
Although these triples will not exhibit the EKM because all objects are equal mass, dynamically formed triples with unequal masses 
produce a similar e Q ct distribution. (See Fig. Il6h 


nents, one of which is approximately thermal for high ei n . 
In the case of triple scattering there is also a strong compo¬ 
nent peaked near ei n ~ 0. This is an artifact of the fact that 
we initialized triples with e 1X] _ = 0 and many new triples were 
formed from an exchange between the tertiary and the inter¬ 
loping star. Since these exchange reactions do not strongly 
affect the inner binary there remains a peak at ei n ~ 0. 

It is clear from this orbital parameter distribution that 
a significant fraction of dynamically formed triples un¬ 
dergo strong KL oscillations. The period of these oscilla- 
tions, £kl, for equal mass systems is given approximately b y 
(H olman et al.|[l997l ;H nnancnet al. I ll997l : lAntogninfll2015l l 

t _ ® ^out (1 _2 \3/ 2 

KL ~ 157T Pin ° U ‘) ' (29) 

Although our calculations are unitless, if we fix the ini¬ 
tial semi-major axes of the binaries to 1 AU in the case 


of binary-binary scattering and in the case of triple scat¬ 
tering we fix the initial semi-major axis of the inner binary 
to 1 AU and set all masses to 1 Mq we can calculate the 
resulting distribution of Ikl in years. This distribution is 
presented in Fig. 1151 There are two interesting features: the 
first is that they are fairly narrow; the vast majority of dy¬ 
namically formed triples span only one decade in fKL. The 
second is that fKL is comparable to the initial dynamical 
timescale of the systems. In the case of binary-binary scat¬ 
tering, the initial periods of the binaries are 0.7 yr and the 
KL timescales are peaked at approximately 3 yr. Likewise 
in the case of triple scattering the outer period is 18 yr and 
the KL timescales have a peak at around 20-30 yr. In the 
case of triple-single scattering there is an additional peak 
at even shorter timescales. Such short KL timescales are a 
consequence of two facts: (1) dynamically formed triples are 
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Figure 15. The distribution of the timescale of KL oscillations of 
dynamically formed triples in a set of model experiments. In the 
case of binary-binary scattering the initial semi-major axes were 
set to 1 AU, and in the case of triple scattering we set a ln = 1 AU 
and a ou t = 10 AU. The masses are all 1 Mq. Triple scattering 
generally produces larger triples and hence longer KL timescales 
because there is a larger scale in the problem, namely a ou t. The 
distributions of 4 kl do not extend much beyond a decade, ex¬ 
cept in the case of triple-binary scattering which has a longer tail 
towards large txL- 


extremely compact, and (2) at the critical velocity triples 
are more easily formed by hardening binaries than by soft¬ 
ening them, thereby leading to smaller semi-major axes and 
shorter periods. 


5.3 Population study 


The above analysis reveals that triples produced from our 
model scattering experiments are extremely compact. It is 
possible that this result is an artifact from the fact that 
in binary-binary scattering we scattered binaries with equal 
semi-major axes and in triple scattering we scattered rel¬ 
atively compact triples. Binary and triple systems in the 
Galaxy span many decades of semi-major axis, so the com¬ 
ponents of typical scattering events will not be of compara¬ 
ble scales. 

To examine the distribution of orbital parameters of 
triples produced from a more realistic population of bi¬ 
nary and triple systems in the Galaxy we run 10 6 apiece 
of binary-binary, triple-single, and triple-binary scattering 
experiments. In each expe riment we draw a p rimary mass 
from the IMF provided bv iMaschbergeij (l2013h which is an 
ana lytic approximat ion to the IMFs found by Kroupal (1200 il l 
and IChabrieij (l2003l l We then draw secondary and, if appli¬ 
cable, tertiary masses from a uniform distribution bounded 
by 0.08 Mq and the primary mass. We next draw semi¬ 
major axes from a distribution uniform in In a from 0.01 to 
10 5 AU. In the case of a triple, we draw two semi-major axes 
and take the smaller to be the semi-major axis of the inner 


binary. We then draw eccentricities from a uniform distr i- 
bution fe.g.. lRaghavan et al.|[2oiol : lDuchene fe KrausH2013h . 
Lastly we determine the stability of the triple by applying 
the Mardling stability criterion (equation 11411 . If the triple is 
unstable we throw out the entire experiment and resample 
all parameters. We set the incoming velocity to 40 km s~\ 
which is approxi mately the velocity dispersion of older stars 
in the thin disk (lEdvardsson et ahlll993l : iDehnen fe Binnevl 

1 19981 : iBinnev fe Merrifieldlll998lb Because we are only inter¬ 
ested in the orbital parameters here, and not the cross sec¬ 
tions for collisions, we assume all stars to be point masses. 

The distribution of orbital parameters of new triples is 
shown in Fig. 1161 Ins tead of e oc t (equa tion 1281) we calculate 
e 0 ct ,m, defined to be jNaoz et al]|2013r ) 


€oct ,M 


{ mm - mn2 
\ mm + TO112 


€oct ? 


(30) 


where mm and mm are the masses of the two components 
of the inner binary. This term accounts for the fact that the 
octupole term in the Hamiltonian has a mass term which 
vanishes in the case of an equal mass inner binary. Thus, 
e oc t ,m indicates the strength of eccentric KL oscillations. 
We estimate the timescale for eccentr ic KL oscillation s (i.e., 
the time between orbital flips) to be (iAntogninill2015l l 


Iekm = -rz - tK t .. (31) 

yj Coct ,M 

These results indicate that dynamically formed triples 
in the field are, as in our model system, extremely compact. 
The peak of the e oc t ,m distribution is broader than the e oc t 
distribution, but centered on ~10~ 2 , indicating that the vast 
majority of dynamically formed triples at high inclination 
undergo strong eccentric KL oscillations. We additionally 
show the KL timescales and eccentric KL mechanism (EKM) 
timescales for these triples in Fig. 1171 Despite the wide range 
of orbital periods in the initial systems (with outer periods 
up to ~10 7 yr), the KL timescales are all quite short. Nearly 
all of them are less than 1000 yr. This is due to two facts: 
(1) the dynamically formed triples are all extremely com¬ 
pact, and so don’t have KL timescales much longer than 
the outer orbital period; and (2) above periods of ~1000 yr 
the binaries become soft and so the cross section for triple 
formation drops dramatically (see Fig. 01. Furthermore, be¬ 
cause e oc t ,m is so large, the EKM timescales are generally 
not substantially larger than the KL time. 


6 DISCUSSION 

Here we apply our results on binary-binary, triple-single 
and triple-binary scattering in four contexts. In Section T6. II 
we compare the production rate of WD-WD binaries with 
highly inclined tertiaries via scattering to the SN la rate. In 
Section IQ we estimate the number of planets in multiple 
systems ejected due to scattering and compare to the esti¬ 
mated number of free-floating planets. We then refine our 
previous calculations by modeling stars with non-zero radii 
to estimate the rate of stellar collisions in Section [60 Fi¬ 
nally, in Section ^.4l we address the question of how long KL 
oscillations persist for triple systems in globular clusters and 
in the field until they are disrupted by scattering events. 
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Figure 16. Orbital parameters of new triples resulting from scattering experiments in our population study (Section l5.3b . The panels 
are as in Fig. Il4l except for the bottom right panel, where we have used e oc t,M instead of e oc t (see equation I30( l. As in the model case, 
triples formed from scattering are extremely compact and have large e oc t- 
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Figure 17. Distribution of the timescales of KL oscillations (left panel) and the eccentric KL mechanism (right panel) in dynamically 
formed triples in our population study (Section l5.3t) . Dynamically formed triples have very short KL timescales, with nearly all KL 
timescales being less than 1000 yr. This is a consequence of the fact that the triples formed are compact and very wide systems are soft, 
and therefore have low interaction cross sections. Because the triples are so compact the timescale for the eccentric KL mechanism is 
also very short, nearly always less than 10 5 yr. 
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6.1 Application to Type la Supernovae 


The progenitor systems of SNe la are currently unknown. 
The two most viable models are the single degenerate 
model, in which a WD accretes mass from a main se - 
quence or giant star (| Whelan fc Iberilll 9731 : iNomo t5 ll982h . 
and the double degenerate model, in w hich two WDs merge 
dlben fe Tutukovlll984l ; IWebbinklil984l 'l. In the past several 
years the double d egenerate model ha s received much obser¬ 
vational support dEU qometjd. 201^_ Schaeferjfc^agnottal 


vational support m joometjaJJ 12LUA; lachactcr & Ragnotta 

201^ Edwardset al] 2012 GonzalezHernandez et al.l20Lfl 


Maoz et al.ll2014l : IShappee et al .11201.'l l, but population syn¬ 
thesis models demonstrate that while it may be possible to 
match the rate of sub- Chandrasekhar mass WD-WD merg- 
ers to the SN la rate (jRuiter et al.ll2009l . 1201 ll : iMaoz et all 
l2012h . hydrodynamical simulations suggest th at a large frac¬ 
tion of these mergers may not produce SNe la (INomoto et alJ 
Il997l : IMaoz fc Mannucc1l2012l : IMaoz et al.ll2014h . 

IThomps'onl (l201lh provided the first analysis of WD- 
WD mergers in triple systems and showed that the per¬ 
turbative influence of the tertiary can decrease the merger 
time by several order s of magnitude llMiller fe Hamiltonl 
120021 : iBlaes et al.ll200 . iThompsonl <1201 il l argued that most 
SNe la may occur in triple systems due to the relatively 
large triple fraction and the larger fraction of triple sys¬ 
tems that merg e rap idl y re l ative to isolated WD-WD bina¬ 
ries. Moreover, IThomps'onl (1201 il l spe culated that s o me o f 
these mergers may result in a collision. iKatz fe Pond (l2012l l 
first showed that clean, head-on collisions of two WDs are 
possible in triple systems using direct N-body calculations 
and claimed that ~5% of WD-WD binaries with tertiaries 
produce a direct collision as a result of non-secular dynam¬ 
ics. The triple scenario for SNe la is attractive because col¬ 
lisions potentially soften the require ment that the merger 
mass exceed the Chandrasekhar mass dRosswog et al 1 120091: 

I Raskin et all l20ld : iKushnir et all 1 201. 'll) . Tertiaries migh t 
also hide WD-WD progenitor systems dThompsonl 1201111 . 
Since the rate of sub-Chandr asekhar mass mergers is 
comparable to the SN la rate (Ivan Kerkwijk et al.l l20ld : 

IMaoz fe MannuccHl2012l : IMaoz et al.ll2012l . l20l4) . this model 

may solve several problems with the double degenerate 
model. 


The triple scenario has an important obstacle, however. 
The same KL mechanism that drives the WD-WD bina¬ 
ries to rapidly merge or collide would drive the stars to 
tidal contact or merger on the main sequence. The main 
sequence stars would then have either tidally circularized 
at a few stellar radii or have undergone common envelope 
evolution. In both cases, the semi-major axis ratio between 
the inner binary and tertiary would likely increase signifi¬ 
cantly, rendering the KL mechanism potentially ineffective 
by the time the stars evolved into WDs. Furthermore, mass 
loss from one of the stars generally exacerbates this problem 
by decreasing the mass ratio of the inner bin ary and trig¬ 
gerin g the eccentric KL mechanism dShappee fe Thompsonl 
l2013lh The eccentric KL mechanis m can drive the inner bi¬ 
nary to much h igher eccentricities dLithwick fe N aoz 2011; 
IKatz et al.ll201ll ); thus, even stars that had escaped tidal cir¬ 
cularization on the main sequence may be brought into tidal 
contact after the primary evolves into a WD. While this may 
be an important source of single degenerate systems, the 
resulting double degenerate system when the smaller mass 


star evolves into a WD may be too wide to merge within a 
Hubble time. It is, however, possible that some fraction of 
these systems might evolve into tight WD-WD binaries dur¬ 
ing common envelope evolution, though it is currently very 
diffic ult to make estimates of this fraction (I Ivanova et al.l 
l2013f) . 

Naively, for the triple scenario to work, it must be pos¬ 
sible for a binary to begin its life in a system that does not 
undergo KL oscillations and then efficiently change into a 
system that undergoes KL oscillations after evolving into a 
WD-WD binary. Scattering is a candidate by which this can 
occur. Although a detailed, comprehensive calculation of the 
rate is left for future work, we can present an estimate to 
determine if scattering is a viable channel for the production 
of high inclination triples. In this estimate we consider three 
environments: the field, open clusters and globular clusters. 
We study the field and open clusters separately because the 
SN la delay time distribution is often decomposed into two 
separate parts: a prompt component which explodes within 
a few hundred Myr of star formation, and a delayed com- 
pone nt which explodes a Gyr or more after star formation 
(e.g., Scamm pjeco&^ildstenl 120051 : iMannucci et all 120061 : 
IMaoz fe Badenesl l20ld h Since the KL timescale is much 
shorter than the scattering timescale (see Section [6.411 . any 
KL-induced collisions may occur rapidly after the scatter¬ 
ing event (e.g., figure 1 of lKatz fe Donell2012l '). The maxi¬ 
mum lifeti me of open clusters is gene r ally only a few hun - 
dred Myr dde La Fuente MarcosHl997l ; lLamers et al.ll2005l h 
so it is plausible that scattering events that contribute to 
the prompt component occur in open clusters and scatter¬ 
ing events that contribute to the delayed component occur 
in the held. However, WD-WD mergers due to KL oscilla¬ 
tions do not necessarily occur instantaneously. Particularly 
in the case of gravitational wave driven mergers, KL os¬ 
cillations can take more than a Gv r to drive the WD-WD 
binary to merger dThompsonll201ll h It is therefore possible 
that scattering events in open clusters could produce SN la 
progenitors which eventually contribute to the delayed com¬ 
ponent of the SN la delay time distribution rather than the 
prompt component. 

We additionally consider scattering in globular clusters. 
The SN la rate in globular clusters is poorly constrained. 
Based on a lack o f d etec tions of SNe la in globular clusters, 

I Voss fe Nelemanjl d2012l) placed an upper limit on the glob¬ 
ular cluster SN la rate, but this limi t was larger than the 
Galactic SN la rate by a factor of 50. iGraham et al.l (|2015l l 
found one hostless SN la which they claimed was likely to be 
in a globular cluster. Based on this detection they estimated 
that the SN la rate might be enhanced in globular clusters 
by a factor of ~25. 

The Galactic SN la rate is Tia ~ 5 x IIP 3 yr” 1 
dMaoz fe Mannuccil 1201 j l. If scattering in the field is the 
primary channel to form high-inclination triple SN la pro¬ 
genitors, then the rate at which scattering events produce 
high-inclination triples must be at least as large as the SN 
la rate. The rate of dynamical formation of new triples may 
be written 


T n ew triple — N TICJV , (32) 

where N is the total number of binaries in the case of binary¬ 
binary scattering or the total number of triples in the case 
of triple scattering, n is the number density of stars, and <r 


© 0000 RAS, MNRAS 000, 000-000 









































































































22 Antognini & Thompson 


is the cross section to produce a new triple. If we consider 
triples or binaries in a small range of semi-major axes, da, 
then the differential rate is 

dV(a) = ^-nva(a) da. (33) 

For our calculation we assume Opik’s law (l6pikHl924T ). 


dN Np 
da a ’ 


(34) 


where No is a normalization constant set by the total number 
of binaries or triples, N: 


N 0 = N 


In 


^max 


Q'min 


-1 


(35) 


Altho ugh a wide log-normal dist ribution is more accurate 
(e.g., iDuquennov fe Mavou ll99lT ). the flat distribution is 
simpler to work with analytically and provides an upper 
limit on the number of systems at the widest semi-major 
axes. 

The cross section for collision can be written in terms 
of d by noting that 

a(a) = na 2 d(v). (36) 


We found in Section EEH that the normalized cross section 
for new triple formation from exchange is proportional to 
v~ 2 for v 1 and is proportional to v~ 6 for v 1. We may 
therefore write the cross section as 


2 - d:ri I 

7T a a 0 —— 


iraaori 
v 2 


and 


v < 1 


(37) 


2" Arit 

na (Jo — n- 
v b 


naorf 

av 6 


v > 1 


(38) 


where do is the normalized cross section for high-inclination 
triple formation at the critical velocity and where we have 
defined rj 2 to be 


V — ^^crit 


(39) 


so as to separate out the dependence of u cr i t on a. rj is the 
same for binaries or triples of any size as long as the masses 
and semi-major axis ratio remain constant (cf. equations 1101 
and E). For the triples we are considering (0.6 Mg and 
a = 10) we have rj 2 ~ 10 3 AU km 2 s~ 2 . 

We additionally found in Section f4.21 that at very high 
velocities the cross section for new triple formation from 
triple-single scattering is dominated by scrambles, which 
have a v~ 2 dependence. This contribution to the overall 
cross section for new triple formation may be written 


2 

2 ~ Aril 

7YCL CT() sc. 7T ~ 
V z 


, naao,sc.'n 


(40) 


where do, S c. is the normalized cross section for scrambles at 
the critical velocity. 

Substituting equations (IM1) . OBI) . (1371) . and (PHI) into 
equation (1331) and integrating with respect to a, we find that 
for binary-binary scattering 


nNondori 


(Ucrit Umin) , 


and 


(41) 


r a > 


nNondori 6 / 1 1 

V 6 \ Ucrit a ma 


(42) 


where we have defined o cr it such that the incoming velocity 
is equal to u cr it, 


Acrit — 


Gmnmi2 

liv 2 


(43) 


For triple-single scattering r a<acrit remains unchanged, but 
we instead have 


nNonrj 2 

v 




T do,sc. Umax , 


(44) 

where we have assumed that a max A a m i n . We may then 
combine these to obtain the overall rate of new triple for¬ 
mation, given by 


r triple — 


it Nona or] 


. f^\ ( 1 1 

(Ucrit U m i n ) T I — 1 I 

/ \ Ucrit Oma: 

for binary-binary scattering, and given by 


(45) 


- new triple 


(^crit ^min) + 


nNondorj 


£1 (-L- 


(d J y a cr it 


< Jo,i 


(Jo 


(46) 


for triple-single scattering. In most cases we have a m in 
a C rit -C a max , so this may be further simplified to 


- new triple 


TtNonaor] 2 


Acrit + 


© 


T]\ 4 1 

Acrit 


(47) 


for binary-binary scattering, and 


r 


new triple 


nNonaor] 2 

V 


X 


Acrit + 



1 

Acrit 



(48) 


for triple-single scattering. This implies that nearly all 
binary-binary scattering events which produce new triples 
occur for systems in which the outer semi-major axis is such 
that the critical velocity of the outer binary is close to the 
velocity dispersion of the surrounding stars. By contrast, 
new triple formation from triple-single scattering events is 
dominated by scrambles in systems near the maximum semi¬ 
major axis. 

Using our scattering experiments from Section r3.2.4l we 
calculate do for binary-binary scattering by calculating the 
cross section for outcomes in which a triple is formed with 
a mutual inclination of 80° < i < 110°. Although the range 
of inclinations that will produce WD-WD collisions is de¬ 
pendent on the semi-major axis (larger systems will pro¬ 
duce collisions for a narrower range of inclinations), the de¬ 
pendence of the required inclination on the semi-major axis 
is not known. We simply assume the in cli nat ion range of 
80° < i < 110° found by iKatz fe Pond (l2012h to produce 
WD-WD collisions to hold for all scales. The cross sections 
we use should therefore be considered to be accurate only at 
the order of magnitude level. We find that for binary-binary 
scattering do = 2 x 1CD 5 . For triple scattering we calcu¬ 
late do in a similar way, but we only select triples which 
were initialized with inclinations less than the Kozai angle 
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( i < 39.2° or i > 141.8°). We find that for triple scatter¬ 
ing (both triple-single and triple-binary) do = 2 x 10 -3 . We 
likewise calculate the cross section for high-inclination new 
triple formation from scrambles at the critical velocity in a 
similar way and find do,sc. = 7x 1CP 5 . 


6.1.1 Estimate for the field 


To estimate F new triple in the field we make the following 
assumptions: (1 ) the total numb er of white dwarfs in the 
Galaxy is 10 10 (lNapiwotzkill2009l ) : (2) the fraction of white 
dwarfs in WD-WD binaries is 10 per cent (lHolberg||2009l) ; 

(3) the fraction of WD-WD binaries with tertiary compan¬ 
ions is the same as the fraction of s tellar binaries with ter - 
tiary companions, i.e., 20 per cent (iRaghavan et al.l l2010h : 

(4) th e volume de nsity of stars in the field is 0.09 Mg 
pc~ 3 dFlvnn et ah 2QQp); (5) the mean star mass is 0.36 
Mg (|Maschbergerl 20131 ): and (6) t he velocity dispersion of 
stars in the thin disk is 40 km s _1 dEdvardsson et aPll993l : 
iBinnev fc Merrifieldlll998l . p. 656). 

Under these assumptions the total number of WD-WD 
binaries is AW = 10 9 , the total number of WD-WD binaries 
with tertiary companions is IVtrip = 2 x 10 8 , and the num¬ 
ber density in the field is n = 0.25 pc -3 . We additionally 
take a max = 10 4 AU and a m i n = 10 -2 AU. From equation 
(|35|) we then have A r o,bin = 7 x 10 7 , and Ao.trip = 10 7 , 
where AWin is the normalization for WD-WD binaries and 
AWrip is the normalization for WD-WD binaries with ter¬ 
tiary companions. Using our data from Section [3.2.41 we cal¬ 
culate the normalized cross section for new triple formation 
to be do = 2 x 10 -3 at the critical velocity of the tertiary for 
triple scattering and do = 2 x 10 -5 for binary-binary scat¬ 
tering. For systems consisting of solar mass stars with an 
incoming velocity of 40 km s _ , a cr it = 0.8 AU. With these 
numbers, equation ([171) gives r new triple ~5x 10~ 9 yr~ 4 for 
triple scattering and F new triple ~ 2x 10“ 13 yr _1 for binary¬ 
binary scattering in the Galaxy. The rates from triple scat¬ 
tering are several orders of magnitude larger than the rates 
from binary-binary scattering because the rate from triple 
scattering is dominated by scrambles at large semi-major 
axes, which are not a possible outcome of binary-binary scat¬ 
tering. However, rates from both binary-binary and triple 
scattering are far below Fi a . In elliptical galaxies the rates 
are further depressed relative to the estimate in the Galactic 
field due to the lower number densities and higher velocity 
dispersions (and hence lower a cr it). 


6.1.2 Estimate for open clusters 

To estimate r new triple for open clusters we make the fol¬ 
lowing assumptions: (1) the typical number density of an 
open cluster is 10 pc“ 3 (IPiskunov et al.ll2007h : (2) the typi¬ 
cal velocity dispersion of an open cluster is 0.3 km s _1 ; (3) 
the binary f raction is 50 per cent and the triple fraction is 
10 per cent (IRaghavan et al.ll2010h : ( 4) the total number o f 
open clusters in the Galaxy is 10 5 (IPiskunov et al.l l2006h : 
and (5) t he typical system is in an open cluster with 300 
members dPorras et al.ll2003l) . 

Under these assumptions the total number of triples in 
open clusters is ~3 x 10 6 and the total number of binaries 
in open clusters is ~2 xl0‘. From equation (1351) we have 


AWin — 10 6 and AWrip = 2 x 10 s . Note that with our as¬ 
sumed velocity dispersion we have a cr it ~ 10 4 AU, which 
is comparable to the semi-major axes of the widest binaries 
observed. Thus in an open cluster environment we are in 
the low-velocity regime even for the widest systems. Since 
scrambles are a dominant channel to produce high inclina¬ 
tion triples only at very high velocities, the contribution of 
scrambles can be neglected in open clusters. 

With these numbers, equation ([471) gives r new triple ~ 
2 x 10 -5 yr~ 4 for triple scattering and T new triple ~ 5x lO -4 
yr _1 for binary scattering in all open clusters in the Galaxy. 
Thus, even in open cluster environments the new triple for¬ 
mation rate is below the SN la rate. These estimates suggest 
that it is difficult for scattering to produce high-inclination 
triples at a rate consistent with the SN la rate. Moreover, 
the estimated rate is dominated by scattering events near 
acrit, which is ~10 4 AU. These separations are only smaller 
than the typical separations between stars in the open clus¬ 
ter by a factor of a few, so our implicit assumption that each 
scattering event proceeds in dynamical isolation is likely un¬ 
warranted at the largest semi-major axes (see Section [6.1.4l) . 


6.1.3 Estimate for globular clusters 


To estimate T new triple for globular clusters we make the fol¬ 
lowing assump tions: ( 1 ) the typical mass of a globular cluster 


is 2 x 10 s M n ( 

Harris 1996T); (2) the typical half-mass radius 

is 3 nc (Harris 

19961): (3) there are 150 globular clusters in 

the Galaxv (Harris 

19961): 14) 0.7 per cent of the systems 


in a g lobular cluster are WD-WD binaries (IShara fc Hurlevl 
120021) : (5) the triple fraction rel ative to the binary frac tion is 
the same as the field, i.e., 0.2 (IRaghavan et al.ll20ld )Fl and 
(6) the typical velocity dispersion of a globular cluster is 6 
km s~ 4 (lHarrisI 1 19961 : IBinnev fc Tremainell2008l . p. 31). 

Under these assumptions the average number density in 
the half-mass radius is ^2500 pc -3 , and the typical separa¬ 
tion between stars is ~10 4 AU. To maintain our assumption 
of dynamical isolation we therefore ta ke q max = 10 3 AU . 
Given a median star mass of 0.36 Mg dMaschbergeil I2OI3I ) . 
the total number of stars within the half-mass radii of all 
globular clusters is ~4x 10 7 , the number of WD-WD binaries 
is ~3xlO s , and the total number of WD-WD binaries with 
tertiaries is ~6x 10 4 . This then implies that lVo,bin. ~ 3 x 10 6 
and lVo,trip. ~ 7 x 10 s . We furthermore have a cr ;t ~ 40 AU. 

With these numbers, equation ([471) gives r new triple ~ 
10' 7 yr~ 4 for binary-binary scattering and r new triple ~ 
4 x 10~ 6 yr~ 4 for triple scattering for all globular clusters 
in the Galaxy. To compare these rates with the Galactic SN 
la rate, we note that the total mass in globular clusters is 
~3xl0' M@ whereas the total s tellar mass in the Galactic 
disk is ~5xl0 10 Mg dMcMillanll201ll) . Thus the fraction of 
stellar mass in globular clusters is ~5x 10~ 4 , so the expected 
SN la rate in globular clusters assuming no enhancement is 


5 In reality the triple fraction in globular clusters will be deter¬ 
mined by the dynamical evolution of the cluster, so there is no 
reason that the triple fraction in globular clusters should be sim¬ 
ilar to that of the field. However, the triple fraction in globular 
clusters is poorly constrained, though it is unlikely to be very 
much higher than that of the field, so we assume here a triple 
fraction to binary fraction of 0.2 as a rough upper limit. 
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~3xlCP 6 yr^ 1 . This is comparable to the scattering rates 
that we find. Thus, if the SN la rate is not enhanced in 
globular clusters, scattering could be an important channel 
to form WD-WD binaries with high inclination tertiaries. 
If, however, t he SN la rate is enha nced by a factor of ~25 
as claimed bv lGraham et all 12015 1 scattering would be un¬ 
likely to be the primary channel by which to produce SN la 
progenitors in triples. 

6 .I .4 Summary of rate estimates 

Although there are substantial uncertainties and simplifica¬ 
tions in the analysis of this subsection, these results imply 
that if the triple scenario is correct, either scattering of very 
wide triples in open clusters leads to the formation of high- 
inclination triple systems or a mechanism other than scatter¬ 
ing is responsible for the formation of high-inclination triple 
systems. It also may be that the full dynamics of open clus¬ 
ters beyond the isolated scattering events we have consid¬ 
ered here (such as the dissociation of the cluster) are impor¬ 
tant for the formation of WD-WD binaries in high inclina¬ 
tion triple systems. While we have treated scattering events 
in open clusters in isolation, this is likely not a good ap¬ 
proximation, particularly at the la rge semi-major axes con¬ 
sidered here (IGeller fc Leighll2015f ). It is also possible that 
triple-binary scattering produces quadruples, which undergo 
strong KL oscillatio ns over a wider range of parameter space 
(IPeicha et al.ll2013f l. Scattering in globular clusters is more 
promising, as the rate of new triple formation can match 
the Galactic SN la rate (after correcting for the fraction of 
mass in globular clusters). Finally, we note that our calcu¬ 
lations of the normalized cross sections were done assuming 
relatively compact triples (semi-major axis ratios of 10) and 
binaries of equal size. The distribution of semi-major axis 
ratios of triples in the Galaxy will certainly be broader than 
we assumed, so the true normalized cross sections will be 
different from the cross sections we calculated. Indeed, they 
will likely be lower due to the fact that scattering is gener¬ 
ally more efficient in compact triples—see Section l5.3l These 
issues should be investigated more fully in a future work. 

We conclude that scattering in the field cannot produce 
WD-WD binaries in high-inclination triples at a rate consis¬ 
tent with the SN la rate, but we cannot rule out the possi¬ 
bility that scattering in open clusters and globular clusters 
leads to high inclination triples at rates consistent with the 
SN la rate. 


6.2 Application to free-floating planets 

The dynamical effect of interloping stars or bin aries on plan- 
etary systems has been studied extens ively (see I Adam 120101 
for a review). iMalmberg et al.l d2007l l measured the rate of 
close encounters between stars in an open cluster environ¬ 
ment and found that dynamical processes will alter the pop- 
ulation of planetary syste ms similar to the Solar System. 
IMalmberg fc Davies! d2009l l found that these close encoun¬ 
ters can produce eccentric orbits with an eccentricity distri¬ 
bution similar to the observ ed distribution of eccen tricities 
of very eccentric exoplanets. iMalmberg et al.l d201l|j calcu¬ 
lated the long-term effect of distant flybys on planetary sys¬ 
tems and showed that ~10 per cent of systems are pushed 


from a configuration which is stable on timescales of ~10 8 
yr to a config uration which is uns table on these timescales. 
Most recently, iLi fc Adamsl (l2015l l calculated the cross sec¬ 
tions for various outcomes and changes to the orbital pa¬ 
rameters from scattering events between a planetary system 
and a binary star. In this subsection we apply our results to 
the free-floating planet population. 

Microlensing surveys have claimed the exis- 
tence of a large populat i on of free - floati ng planets 
dZapatero Osorio et aid 120001 : ISumi et akl l201ll l . It is un¬ 
known whether these sub-brown dwarf objects form in 
isolation or if they were ejected from the planetary systems 
in which they were born. Numerical simulations indicate 
that it may be difficult for planet-planet scattering to fully 
account for the ob s erved numbers of free-floating planets 
dVeras fc Ravmondl l2012l l. While there is some observa¬ 
tional evidence for the isolated for mation of su b-brown 
dwarf objects in the Rosette Nebula dGahm et al.l l2013| j. it 
is nevertheless unclear if isolated formation is the dominant 
channel for the production of free-floating planets. 

We here examine the efficiency of planet ionization due 
to scattering events. We perform three classes of scattering 
events: (1) a multi-planet system scattering off of a single 
or binary star system, (2) a circumbinary planetary system 
(i.e., a P-type orbit) scattering off of a single or binary star 
system, and (3) a planetary system in a binary (i.e., an S- 
type orbit) scattering off of a single or binary star system. 
We furthermore examine planet scattering in both open clus¬ 
ters and in the field. In these experiments the masses of the 
stars are set to 1 Mg and the planets are set to 0.01 Mg. 
In the triple systems the inner object is given a semi-major 
axis of 1 AU and the outer object is given a semi-major axis 
of 20 AU and the orbits are set to be coplanar. (In the case 
of a multi-planet system we instead set the semi-major axis 
of the inner planet to 4 AU to encourage planet-planet scat¬ 
tering.) In the case of scattering off of a binary, the incoming 
binary is given a semi-major axis of 100 AU. For the field 
the incoming velocity is set to 40 km s -1 and for open clus¬ 
ters it is set to 3 km s -1 . Our assumed incoming velocity 
for open clusters is larger in this section than in Sections l6.ll 
and IOI due to the fact that the scattering experiments we 
perform in this subsection are prohibitively expensive at an 
incoming velocity of 0.3 km s . Consequently we calculate 
the cross sections at this larger incoming velocity and later 
scale the results to a lower incoming velocity. 

With these initial conditions we then calculate the cross 
section for planet ionization and present them in Table El 
In all systems in the field the normalized cross section is 
within an order of magnitude of unity. In open clusters the 
cross section is enhanced by one to two orders of magni¬ 
tude because the incoming velocity is smaller by a factor of 
~10. Multi-planet-single scattering has the largest normal¬ 
ized cross section due to induced planet-planet scattering. 
Although this effect is also present in multi-planet-binary 
scattering, its normalized cross section is much lower simply 
due to the much larger area of the interloping binary. The 
physical cross section for planet ionization is slightly larger 
in multi-planet-binary scattering than multi-planet-single 
scattering. 
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Table 6. Cross sections for planet ionization. We include both 
the normalized cross section, <x, and the physical cross section, a. 
Multi-planet systems consist of two planets, circumbinary systems 
consist of a planet orbiting a binary, and S-type planets consist 
of a planet orbiting a star in a binary system. 


Field environment 


System 

a 

<7 (AU 2 ) 

Multi-planet—single 

6.52 ± 1.14 

2608 ± 457 

Multi-planet-binary 

0.32 ± 0.05 

3325 ± 508 

Circumbinary-single 

3.21 ± 0.08 

1285 ± 31 

Circumbinary-binary 

4.00 ± 0.09 

41648 ± 950 

S-type planet-single 

0.68 ± 0.02 

272 ± 7 

S-type planet-binary 

1.26 ± 0.03 

13147 ± 309 

Open cluster environment 

System 

<T 

cr (AU 2 ) 

Multi-planet—single 

68.73 ± 14.57 

27493 ± 5828 

Multi-planet-binary 

15.57 ± 4.44 

161879 ± 46219 

Circumbinary-single 

53.59 ± 6.69 

21436 ± 2675 

Circumbinary-binary 

15.20 ± 3.51 

158082 ± 36466 

S-type planet-single 

64.43 ± 2.03 

25771 ± 811 

S-type planet-binary 

7.57 ± 0.77 

78703 ± 8046 


The timescale for planet ionization can be written 


^scatter — 


1.1 x 10 

a 


1 

nnaal ut v 



-l 


/ Qout \~ 2 / V 

\20AU/ \3kms -1 


yr- 

(49) 


Even for the outcomes with the largest cross sections (<r ~ 
70) the ionization timescale for a given system in a cluster 
is ~20 Gyr. We noted above that a was calculated for an 
incoming velocity of 3 km s _1 . For a more realistic incom¬ 
ing velocity of 0.3 km s“ , a would increase by two orders of 
magnitude due to gravitational focusing. However, the scat¬ 
tering timescale is also proportional to u _1 , so the overall 
effect of a smaller incoming velocity would be to decrease the 
ionization timescale by one order of magnitude to ~2 Gyr. 
This implies that planet ionization from scattering with in¬ 
terloping stellar systems is rare in stellar clusters; only ~10 
per cent of systems with planets on wide orbits would have 
lost them in a cluster with an age of 200 Myr. In the held, 
planet ionization occurs even less frequently due to the lower 
densities and cross sections. Assuming a stellar density of 
0.1 pc -3 , a typical velocity of 40 km s -1 , and a normalized 
cross section of ~6, the scattering timescale is ~4x 10 12 yr. 
Thus, fewer than 1 per cent of systems with planets on wide 
orbits would have lost them in the lifetime of the Galaxy. 
Scattering from interloping systems can therefore produce 
some free-floating planets, but not at the level claimed by 
ISumi et al.l d201ll ). 


6.3 Stellar collisions during scattering events 

Although the distances between stars are typically much 
larger than their radii, close encounters between stars can 
result in tidal interactions, mergers, and collisions. Can¬ 
didates for observed outbursts from collisions or merg¬ 


ers include V838 Mon and V1309 S co (iBo nd et al. 2003) 
Tvlenda fe Sokcr l l2006l : lTvlenda et al.ll2011 : Kochanek et al.l 


2014). It is still unclear whether they should prefer¬ 
entially occur between stars born in the same system 
or between stars in d ifferent systems dLeigh et al.l l201ll : 
IPerets fe KratteH 120121 ). If collisions are typically between 
stars in the same system, then these collisions would be 
driven either by dynamical instabilities in the primordial 
system (in which case collisions will occur not long after 
the protostar reaches the main sequence, if it ever does) or 
due to KL oscillations. If, however, a substantial fraction of 
collisions are between stars in different systems, then this 
necessarily implies that scattering would be an important 
channel for stellar collisions. 


6.3.1 Numerical calculations 

Collisions are always a possibility in resonant scattering 
events because such events are chaotic. We here rerun our 
model system (ai n = 1 AU, a ou t = 10 AU—see Section GEE) 
but we now give the stars radii of 1 Rq. We do not include 
the effects of tides. We halt the calculation after a single 
collision, so multiple collisions are not considered. 

To calculate these cross sections we must account for 
the fact that some fraction of the observed collisions in 
triple-single and triple-binary scattering are due to KL os¬ 
cillations and are not induced by scattering. To correct for 
this effect we run an identical set of triples in isolation for 
the same length of time and calculate the fraction of triples 
that collide. This fraction (~few per cent) is then subtracted 
from the observed collision fraction from the scattering ex¬ 
periments before the cross section is calculated. We find 
that the normalized cross section for collision at the crit¬ 
ical velocity is do = 0.109 ± 0.004 for binary-binary scat¬ 
tering, do = 0.134 ± 0.004 for triple-single scattering, and 
do = 0.235 ± 0.007 for triple-binary scattering. 

Because these calculations are all done in the equal mass 
case, they are biased towards scattering-induced collisions 
and away from KL-induced collisions because EKM oscilla¬ 
tions will be stronger, leading to a larger collision fraction 
for triples in isolation. It is also plausible that at more ex¬ 
treme mass ratios ionizations will be a more common out¬ 
come (see Section r3.2.3ll . leading to fewer scattering-induced 
collisions. This implies that the collision cross sections de¬ 
rived above are likely upper limits. We do not calculate the 
velocity dependence of the collision cross sections, however, 
so we cannot rigorously assert that these cross sections are 
upper limits. 


6.3.2 Rate calculation 

We show in Appendix[A]that the velocity dependence of the 
collision cross section for three-body interactions scales as 
v~ 2 for i < 1 and v~ 6 for v 1 so long as v is much less 
than the escape speed of the stars. Collisions are therefore 
an exchange-like process. Because of this, the rate of stel¬ 
lar collisions from scattering may be estimated in the same 
way as the new triple formation rate was estimated in Sec¬ 
tion RH1 We here follow the same analysis to estimate the 
rate of collisions due to scattering with one modification. Be¬ 
cause we hold the stellar radii fixed, the ratio ( R/a) varies 
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as we vary a. iFregeau et al.l d2004l') found that the collision 
cross section scales linearly with the ratio ( R/a ) for binary¬ 
binary scattering, so we incorporate this factor into the rate 
calculation: 


and 


2 ^ h:ri 1, 

u = 7ra (Jo — 7t 

v 


' a o\ _ napaor] 


a) 


v < 1, 


(50) 


2- ^crit /O 0\ 

1 7ra a ° - 6 ~ l - ) — 

v° \ a / 


naocroV 
a 2 v 6 


.6 


v > 1, 


(51) 


where ao is the outer semi-major axis at which do was calcu¬ 
lated of 10 AU and rj is defined as in Section 16.11 Note that 
we assume here that the collision cross sections for triple¬ 
single and triple-binary scattering carry the same (R/a) de¬ 
pendence as binary-binary scattering. Proceeding to inte¬ 
grate equation (EH) as we did before, we find 


r collision — 


7rNonaor] 2 ao 


In 


f &crit 
\ &min 


+ 


1 / 7/ 


(52) 

where, as in equation (S3. we have neglected a max . 

We now evaluate r co iiision for the field, open cluster 
environments, and globular cluster environments using the 
same assumptions as we did in Section r6.ll We additionally 
assume that the mass o f the Galactic disk is 5 x 10 10 Mq 
(Dclinen & Binn cv 1998|). Given a mean stellar mass of 0.36 


Mq (lMaschbergeiil2013f) . this implies that the total number 
of systems in the Galactic disk is 1.4 x 10 11 . Making use 
of our assumption tha t 10 per cent of systems are triples 
(IRaghavan et al.ll2010r i. this implies that the total number 
of triples is 1.4 x 10 10 , and the normalization constant of 
equation (13511 is No = 10 9 for the field. We assume the same 
normalization constant for open clusters as we did in Sec¬ 
tion IQ of No = 2 x 10 5 . Furthermore, assuming a typical 
stellar mass of 0.36 Mq we find rj 2 ~ 850 AU km 9 s -2 . 
With these numbers we estimate F co ni s i on ~ 2 x 10~ 6 yr -1 
in the field, ~9xl0 -6 yr -1 in open cluster environments. 
For globula r clusters we assum e that the binary fraction is 
3 per cent (iMilone et al.l [2012l l and the ratio of the triple 
fraction to the binary fraction is 0.2 as it is in the field. 
This implies a normalization constant for globular clusters 
of No = 3 x 10 6 and a collision rate of ~8xl0 -4 yr -1 in 
globular cluster environments. The ratio between the esti¬ 
mated rate in the open clusters and the field is much smaller 
here than in the case of our SN la rate estimate because o cr it 
is much larger in open clusters, leading to smaller values of 
(R/a), and hence smaller collision cross sections. 

These rates are much smaller than the observed rate of 
stella r mergers in the Galaxy of ~0.5 yr -1 dKochanek et al.l 
l2014l l. This implies that scattering is a relatively minor con¬ 
tributor to the stellar merger rate. Of course, since ~10 per 
cent of all triples are formed at inclinations |cosi| ^ 0.1, 
many of these systems will undergo merger events, and will 
therefore likely be a much more important contributor to 
the stellar merger rate. We also note that as in Section [6.II 
our assumption of dynamical isolation is likely invalid in an 
open cluster environment. We save a treatment of the full 
dynamics of open clusters for a future work. 


6.4 How long do high inclination triples survive? 

Although KL oscillations can drive interesting behavior in a 
diverse set of systems in isolation, real systems are not iso¬ 
lated. Perturbations to the gravitational potential of a hier¬ 
archical triple generally suppress the KL resonance. While 
several of these effects have already been studied in detail 
(e.g., relativistic precession), scattering has not. 

Scattering will suppress KL oscillations on the timescale 
of an interaction in which the triple is scattered from a high- 
inclination state to a low-inclination state or some other con¬ 
figuration in which KL oscillations are not present (e.g., two 
binaries). Moreover, if this timescale is much shorter than 
the timescale of KL oscillations themselves, KL oscillations 
will not occur at all. 

To calculate the lifetime of high inclination triples we 
performed a suite of scattering experiments over a range of 
outer semi-major axes with a semi-major axis ratio of 10. 
We performed experiments in two conditions: the field, for 
which we take the incoming velocity to be 40 km s _ ; and 
globular clusters, for which we take the incoming velocity 
to be 6 km s . The scattering timescale is calculated from 
the cross section for the triple to either disrupt or change its 
inclination from the range 80° < i < 110° (see Section [6.Ill 
to less than the Kozai angle (i < 39° or i > 141°). We take 
the number density of stars to be 0.25 pc~ 3 in the field and 
10 4 pc~ 3 in globular clusters. 

The scattering timescales for triple-single and triple¬ 
binary scattering, along with fxL is shown in Fig. 1181 The 
scattering timescale, t s cat, is proportio nal to a - 1 because, for 
fixed incoming velocity, a/a 2 oc v 2 lit (iHut fe Bahcalllll983ll 
and v 2 rit oc a. 

In the field the density of stars is so low that high in¬ 
clination triples persist for a Hubble time except for the 
very widest triples (a ou t > 10 5 AU). Globular clusters are 
dense enough, however, that high-inclination triples wider 
than ~few AU will be disrupted in a Hubble time. But 
even in globular clusters, triples as wide as 1000 AU will 
undergo many KL oscillations before disruption by scatter¬ 
ing. This analysis ignores the cumulative effect of many per¬ 
turbations to the angular momentum of the outer binary, 
however. At velocities much larger than the critical veloc¬ 
ity, perturbations to the angular momentum of a binary be¬ 
come much larger than perturbations to its energy. Very 
distant passages cause the angular momentum of wide bi¬ 
naries to undergo a random walk and can generally lead 
to extremely large ecce ntricities before disruption occurs 
dKaib fe Ravmondll2014l l. In the case of a triple with a wide 
outer binary, these perturbations can drive the triple to in¬ 
stability in less time than it would take the triple to disrupt 
due to scattering alone. Moreover, if the cumulative effect 
of many distant passages can be modeled as a global tidal 
field, the outer binary of a wide triple will u ndergo eccen¬ 
tricit y oscillations similar to KL oscillations dKatz fc Dong 
Uml). Since the cumulative effect of multiple distant per¬ 
turbations can have a substantial impact on the orbital pa¬ 
rameters of a triple system on long timescales, particularly 
in dense environments like globular clusters, these long-term 
effects should be more fully investigated. 
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Figure 18. The lifetime of high-inclination triple systems in glob¬ 
ular clusters and in the field. Triple systems in the field persist 
for a Hubble time out to widths of 10 5 AU and triple systems in 
globular clusters persist for a Hubble time (black dotted line) out 
to widths of a few AU. We also include the KL timescale (black 
dashed line). For all but the widest triples in globular clusters the 
scattering time is much longer than £kl- All but the very widest 
triples in globular clusters therefore undergo many KL oscillations 
before disruption by scattering. The scattering timescale is pro¬ 
portional to the inverse of the outer semi-major axis (blue dotted 
line). 


7 CONCLUSIONS 

We have explored the properties of binary-binary, triple¬ 
single, and triple-binary scattering events in a wide vari¬ 
ety of contexts using over 400 million numerical scattering 
experiments. We have calculated the cross sections for the 
outcomes of these scattering events in several model systems 
and grouped these outcomes into a small number of out¬ 
come classes (Tables |T] [7] EH and [3). These outcome classes 
can, in turn, be broadly grouped into “exchange-like” and 
“ionization-like” outcomes based on their velocity depen¬ 
dence in analogy to binary-single scattering. Ionization-like 
outcomes exhibit a v~ 2 dependence at large v and exchange¬ 
like outcomes exhibit a steeper velocity dependence at large 
v (often v~ 6 , but sometimes v ~ 4 —see Fig.jGj. The statistical 
uncertainties on these cross sections are on the order of a few 
per cent or less. The systematic uncertainty is comparable, 
but only serves to increase the cross sections of ionization¬ 
like outcomes. 

We find that the cross section for new triple formation is 
“exchange-like,” and at the critical velocity the normalized 
cross section for new triple formation is do = 0.1 for both 
triple-single and triple-binary scattering, and a o = 0.02 for 
binary-binary scattering. At high velocities the new triple 
cross section has a v~ 6 dependence for binary-binary scatter¬ 
ing. For triple-single scattering, the new triple cross section 
also has a v~ 6 dependence for velocities slightly in excess of 
the critical velocity, but at high velocities, scrambles, which 
have a v~ 2 dependence, are the dominant channel to form 


new triples. Similarly, the cross section for quadruple forma¬ 
tion from triple-binary scattering is an exchange-like process 
with a v~ & dependence at high velocity. These results imply 
that quadruple formation and triple formation in environ¬ 
ments dominated by binary-binary scattering are efficient 
only if the velocity is below the critical velocity (i.e., in 
cluster environments). In high velocity environments, triple¬ 
single scattering is the dominant channel to produce new 
triples. 

We provide analytic fits for the velocity dependence of 
the ionization cross sections in binary-single, binary-binary, 
and triple-single scattering (equations 1201 and 1221) . We also 
provide analytic fits for the velocity dependence of the new 
triple cross sections in binary-binary and triple-single scat¬ 
tering (equations [25] and 1261 . 

We measure the dependences of the cross sections on 
various orbital parameters. We find that the dependence of 
the cross sections on the semi-major axis ratio can be at¬ 
tributed entirely to changes in the critical velocities of the 
various components, at least as long as the semi-major axis 
ratio is greater than ~10 (Fig. [3]. 

We find that there is no eccentricity dependence on the 
cross sections, except insofar as triples with more eccentric 
tertiaries come closer to the boundary of stability (Fig. []]). 
Weaker perturbations to the eccentricity of the tertiary can 
therefore cause the triple to destabilize, leading to a modest 
increase in the ionization cross section for large outer eccen¬ 
tricities. This increase is well modeled using our fit to the 
cross section for changes to the outer eccentricity in flyby 
outcomes. We also studied the mass dependence of the cross 
sections and find that in general cross sections increase for 
ionization-like outcomes in binary-binary and triple-binary 
scattering as the mass of a single star tends toward infin¬ 
ity (Fig. [5}. In triple-single scattering the cross sections for 
all outcomes decrease in the high mass limit because the 
interloping system cannot be disrupted. 

Almost all scattering events are simple flybys which 
leave the hierarchical structure intact. We calculate the cross 
section for changes to the orbital parameters after flybys in 
triple-single and triple-binary scattering events (Fig. HOI) . 
We find that changes to the eccentricity of the tertiary are 
weil fit by a Gompertz function (equation 1271) . The normai- 
ized cross section for a iarge change in the inclination of the 
triple (Acosi ~ 0.5) is <r ~ 0.1 after correcting for changes 
to the inclination due to Kozai-Lidov oscillations (panel c of 

Fig. [TUI. 

We study the properties of triples formed from scat¬ 
tering events and find that dynamically formed triples are 
extremely compact (Figs. [Til and [TUI . Indeed, most dynam¬ 
ically formed triples are very close to the stability bound¬ 
ary although longer-term simulations indicate that they are 
stable for at least hundreds of orbits. Because these triples 
are so compact, the timescale for Kozai-Lidov oscillations 
is very short, sometimes just a few times the outer orbital 
period (Figs. [T5l and [T7I . These triples should therefore ex¬ 
hibit strong non-secular dynamics. We furthermore find that 
the inclination distribution of dynamically formed triples is 
approximately uniform in cos i and the distribution of outer 
eccentricities is approximately thermal. These results hold 
both for a model system and for a population study where 
we model scattering events in the field. 

Because many-body gravitational dynamics are ubiqui- 
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tous in astrophysics, our results have implications in a vari¬ 
ety of different contexts. We find that scattering in the field 
cannot produce high inclination triple systems at a rate con¬ 
sistent with the SN la rate, particularly in ellipticals. This 
result poses problems for the triple scenario for SNe la be¬ 
cause high inclination triples need to be produced after the 
stars that eventually comprise the inner binary of the triple 
evolve into WDs. If this does not occur KL oscillations will 
drive the inner binary to tidal circularization and will sup¬ 
press future KL oscillations. We show in Section eh that 
scattering in open clusters is more efficient and can produce 
high inclination triples at a rate a factor of ~20 below the 
SN la rate. However, it is unlikely that our assumption of 
dynamical isolation in open clusters is valid, so the forma¬ 
tion of high inclination triples may be more efficient than 
we have estimated when the full dynamics of the cluster are 
considered. A triple scenario where intermediate mass stars 
in binaries or triple systems evolve to WDs and then scat¬ 
tering leads to high inclination tertiaries may therefore be 
an important contributor to the prompt component of the 
SN la delay time distribution. If the time for KL oscillations 
to drive the WD-WD binary to merger is much longer than 
the lifetime of the cluster, scattering in open clusters could 
also contribute to the delayed component of the delay time 
distribution. We save a complete exploration of scattering 
in cluster environments and its implications for SNe la for 
a future work. 

Scattering can also lead to collisions between stars. We 
show that the velocity dependence for collisions in binary- 
single scattering due to three-body interactions is exchange¬ 
like with a high-velocity dependence of fi -6 . We estimate 
the rate of stellar collisions in the Galaxy due to triple 
scattering to be ~2xlCP 6 yr -1 . In open cluster environ¬ 
ments the estimated collision rate is comparable, ~9xlCP 6 
yr -1 . In globular cluster environments the estimated colli¬ 
sion rate is larger, ~8xl0 -4 yr -1 , due principally to the 
higher number densities and higher velocity dispersion rel¬ 
ative to open clusters. However, our assumptions about dy¬ 
namical isolntion_arelikeIy not valid in an cluster environ¬ 
ments (iGeller fe Leighll2015l ). The collision rate in clusters 
may therefore be much higher than our naive analysis sug¬ 
gests. 

We additionally apply our results to several different 
types of planetary systems and find that scattering from 
external encounters is a negligible contributor to the free- 
floating planet population compared to planet-planet scat¬ 
tering. However, in dense open clusters nearly 10 per cent 
of planets on wide (~100 AU) orbits will be ionized due to 
scattering with interloping systems. 

We finally examine the stability of KL oscillations to 
scattering events (Fig. 1181) . We find that the cross section 
for the triple to move from a high-inclination to a low- 
inclination state is small enough that the timescale to scatter 
out of a regime in which KL oscillations take place is much 
longer than the Hubble time for systems in the field and 
all systems with outer semi-major axes <few AU in globu¬ 
lar clusters. Furthermore, in both the field and in globular 
clusters the KL timescale is much shorter than the scatter¬ 
ing timescale. Thus triples with moderate semi-major axis 
ratios in any environment may be considered sufficiently iso¬ 
lated that at least several KL oscillations will proceed undis¬ 
turbed. 
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APPENDIX A: THE COLLISION CROSS 
SECTION 

The dependence of the collision cross section in few-body 
scattering problems on the incoming velocity and stellar ra- 
dius has been inve s tigated by a numbe r of studies (e.g., 
iFreeeau et al.l 120041 ; iLeigh fe Gelled 120151 ') , but an analytic 


fit to these dependences does not exist in the literature. We 
present such a fit here in the case of equal mass stars. 

Collisions can occur in two different ways: through a 
three-body interaction or through a two-body interaction. 
In a three-body interaction the interloping star exchanges 
momentum with one star of the binary before eventually 
colliding with the other. In a two-body interaction the in¬ 
terloping star collides into one of the stars of the binary 
without interacting with the other at all. 

For three-body interactions, the collision cross section 
scales with velocity in the same way that the cross section 
for exchange does. That is, for v < v CT it the velocity scal¬ 
ing is v~ 2 and for v > v C rit the velocity scaling is v~ 6 . The 
low-velocity scaling is due to gravitational focusing, whereas 
the high-velocity scaling is likely because the mechanism by 
which collisions occur in three-body scattering events is sim¬ 
ilar to the mechanism by which exchanges occur. Exchanges 
occur whe n two stars undergo a close approach and exchange 
momenta (jHutlll983h . This ejects one star formerly in the 
binary and replaces it with the interloping star on a simi¬ 
lar orbit. A collision may occur if the momentum exchange 
places the interloping star on an orbit eccentric enough to 
lead to a collision. Since the mechanism for collision is sim¬ 
ilar, we would then expect the collision cross section to dis¬ 
play the same v~ 6 velocity dependence as the exchange cross 
section. Note, however, that this is simply a plausibility ar¬ 
gument. We do not attempt here to show rigorously that this 
is the mechanism responsible for the observed velocity de¬ 
pendence. Nevertheless, the numerical calculations indicate 
that the cross sections indeed show a 0 -6 velocity depen¬ 
dence in this regime just like exchange. 

We may combine these velocity dependences as we did 
in Section rm to find the overall cross section: 


1 

O 3-body 



(Al) 


where a a and dt are constants of order unity that must be 
fit to numerical data. 

Let us consider now two-body interactions. Like three- 
body interactions, these exhibit two velocity regimes. How¬ 
ever, in this case the two velocity regimes are separated by 
the escape velocity of the stars, v eB c■ For v v eS c, gravi¬ 
tational focusing from the individual stars is negligible and 
the cross section for collision will approach a constant value 
equal to the geometric cross section of the stars. For veloc¬ 
ities less than v esc , gravitational focusing from individual 
stars will lead to a v~ 2 dependence. For two-body interac¬ 
tions, then, the overall cross section will be 


~ _ * . _L L 

<3"2-body — ^geometric i To •> 


where the geometric cross section is just given by 

-2 

^geometric "7T l 


(A2) 

(A3) 


and a c must be fit to numerical data, though it will be of 

Order G geometricU esc . 

The collision cross sections for two- and three-body in¬ 
teractions may then be combined additively to yield the 
overall collision cross section: 


deoil. — ^2-body T d,3_body — 
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Figure Al. The cross section for a collision to occur in binary- 
single scattering as a function of the normalized incoming velocity. 
We show both numerical data (gray points) and a fit to equation 
JA4H (black dashed line). This system consists of three 1 Mq 
stars each with a radius of 0.3 Rq. The binary orbit is circular 
with a semi-major axis of 1 AU. The best-fitting parameters are 
a a ~ 0.475, ~ 2.00 x 10 -2 , <r c ~ 2.25 X 10 -2 , and ^geometric ~ 

1.61 x 10 -5 . 


A fit of equation (lA4l) to numerical data is shown in Fig. EH 
There is an excellent match between the numerical data and 
the analytic fit. 
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